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I.  INFORMATION  SYSTEMS 


A.  Project  6151.  STATISTICAL  DATA  PROCESSING  AND  PATTERN  RECOGNITION 

Principal  Investigator:  T.  M.  Cover 
Staff:  R.  King,  A.  El  Gamal,  K.  Jarett, 

J.  van  Campenhout , M.  Salehi 

1 . Objective 

The  purpose  of  this  project  is  to  Investigate  problems  of  In- 
formation processing  and  pattern  recognition. 

2 . Current  Status  of  Work 

a . The  Measurement-Selection  Problem  In  Regression  Analysis 

(J.  van  Campenhout) 

The  similarity  between  measurement  selection  in  linear 
regression  and  In  two-class  common  covariance  gausslan  classification 
problems  has  been  developed.  Based  on  earlier  results  obtained  from 
measurement  selection  In  classlf Ication  problems,  it  can  be  concluded 
that  optimal  subset-selection  algorithms  must  be  exhaustive  in  regres- 
sion problems  with  vector-valued  regressors. 

It  is  also  proved  that  the  subclass  of  regression  problems 
with  scalar  regressors  cannot  Induce  all  monotone  orderings  on  regressor 
subsets.  It  Is  possible,  therefore,  that  none:diaustlve  optimal  search 
algorithms  exist  for  this  class  of  problems. 

b.  Resolution  of  the  Hughes  Paradox  (J.  van  Campenhout) 

Various  researchers  have  observed  that.  In  some  cases,  the 
Bayes  recognition  accuracy  of  pattern  classifiers  peaks  with  Increasing 
measurement  complexity.  Tills  result  Is  apparently  paradoxical  because 
more  complex  measurements  should  reduce  the  mlsclasslf Ication  probabil- 
ity. 

This  problem  has  been  Investigated , and  It  has  been  veri- 
fied that  the  paradoxical  peaking  effect  Is  caused  by  Incorrect  compar- 
isons of  statistically  Incomparable  models.  A formalization  of  compar- 
ability has  been  Introduced,  some  of  the  results  obtained  In  the  liter- 
ature have  been  analyzed,  and  a paper  will  be  published. 
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c.  Practical  Application  of  Superposition  Codes  (K.  Jarett) 

Preliminary  results  indicate  that  superposition  codes  re- 
quire roughly  the  same  decoding  complexity  as  do  random  codes,  given  a 
fixed  probability  of  error  requirement. 

Superposition  codes  are  attractive  because  their  complex- 
ity grows  only  linearly  with  rate.  Instead  of  using  a rate  R random 

code  with  2”^  codewords,  an  L-stage  superposition  code  in  which  L = 

, nRn 

R/Rq  can  be  applied,  and  each  stage  has  R^  and  2 ^ codewords.  The 

decoder  need  only  be  capable  of  L consecutive  maximum-likelihood  de- 

nRf, 

codings  in  a codebook  of  2 The  problem  is  that,  for  a given  block 

length  n,  the  superposition  code  has  a much  higher  error  probability. 

For  large  n,  exponential  error-probability  bounds  in 
the  form  of  2 ” ^ show  that  a 1-stage  R = ^ R^  superposition  code 
requires  a block  length  of  nL  to  equal  the  error  probability  of  a block 
length  n rate  R random  code.  For  this  superposition  code, 

e - ^ 

if  = C/L  and  R^  = R/L  for  all  1;  otherwise,  the  error  bound  is 
worse.  With  a block  length  nL  and  Rj^  = R/L,  each  stage  of  the  super- 
position code  has  a codebook  of  = 2*^^,  As  a result,  for  either 

superposition  or  ordinary  random  codes  (given  a fixed  Pg) > decoder 

must  be  capable  of  a maximum-likelihood  decoding  in  a codebook  of  2^^. 

Numerical  results  obtained  for  the  binary  symmetric  chan- 
nel confirm  these  theoretical  estimates. 

d . Multiple-Access  Channels  with  Generalized  Feedback 
(R.  King) 

An  achievable  rate  region  has  been  determined  for  multi- 
ple-access channels  with  generalized  feedback — feedback  not  necessarily 
equal  to  but  correlated  with  the  channel  output.  Noiseless  feedback  in- 
creases the  capacity  of  a multiple-access  channel  and  is  a special  case 
of  generalized  feedback. 
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The  theory  of  coding  with  side  Information  was  used  to 
extend  noiseless  feedback  to  generalized  feedback.  The  new  achievable 
rate  region  yields  the  total  cooperative  rate  region  for  several  non- 
trivial examples. 

e.  Multiple-Access  Channels  with  Correlated  Sources  and 

Feedback  (R.  King) 

Consider  a multiple-access  channel  with  sources  correlated 
In  a special  manner.  Each  transmitter  observes  one  Independent  source 
(as  In  the  standard  multiple-access  channel);  In  addition,  however,  both 
transmitters  observe  a third  source  In  common,  and  the  capacity  of  this 
channel  Is  a known  result . 

An  achievable  rate  region  has  been  derived  for  this  chan- 
nel with  noiseless  feedback  added.  New  Intuition  Into  the  relationship 
betwwsen  the  known  results  for  multiple-access  channels  with  correlated 
sources  and  with  feedback  Is  thus  gained. 


f . Bounds  on  Cardinality  of  Auxiliary  Random  Variables  In 
Multiuser  Communication  Channels  (M.  Salehl) 


The  method  of  Ahlswede  and  Komer  has  been  employed  to 
determine  the  upper  bounds  on  the  cardinality  of  auxiliary  random  vari- 
ables Involved  In  the  description  of  capacity  regions  for  multiuser  chan- 
nels. Some  of  the  existing  bounds  have  been  Improved. 

Bounds  have  been  successfully  obtained  In  the  following 


cases : 


• source  coding  with  side  Information 

• more  capable  broadcast  channels 

• broadcast  channels  with  degraded  message  sets 

• multiple -access  channels  with  correlated  sources 

• multiple-access  channels  with  feedback 

• common  Information  of  two  random  variables 
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g . Capacity  of  the  Product  and  Sum  of  Two • Inconsistently 
Degraded  Broadcast  Channels  (A.  El  Gamal) 

Let  p(y^,Zj^|x^)  = pCy^  |xj^)  p(z^  |y^)  and  p(y^,z^\x^)  = 
p(z  |x  ) p(y  |z  ) be  two  degraded  broadcast  channels.  The  broadcast 
channel  with  sender  receivers  and  ^i>^2  referred 

to  as  the  product  of  two  Inconsistently  degraded  broadcast  channels. 

Similarly,  the  sum  of  two  Inconsistently  degraded  broad- 
cast channels  Is  defined  as  the  broadcast  channel  with  U Xg  as  sender 
and  Y,  U Y_  and  Z u Z as  receivers . In  contrast  to  the  two  degraded 

1a  1a 

broadcast  channels , the  sum  and  product  channels  are  not  degraded . 

Capacity  regions  have  been  established  for  the 

• discrete  memoryless  product  of  two  Inconsistently  degraded 
broadcast  channels 

• spectral  gausslan  broadcast  channel 

• discrete  memoryless  sum  of  two  Inconsistently  degraded 
broadcast  channels 


These  capacity  theorems  Include  the  partial  results  obtained  by  Poltyrev 
for  the  discrete  memoryless  product  channel  and  by  Hughes -Hartogs  for  the 
spectral  gausslan  broadcast  channel,  and  they  prove  that  the  Hughes-Har- 
togs  rate  region  Is  optimal  for  the  zero  common-rate  case.  A paper  has 
been  submitted  to  Problemy  Peredacl  Informaccl. 


h.  Capacity  of  a Class  of  Broadcast  Channels  (A.  El  Gamal) 

The  capacity  region  has  been  determined  for  any  discrete 
memoryless  broadcast  channel  p(y,z (x)  for  which  I(X;Y)  > I(X;Z)  for 
all  Input  distributions.  This  region  resembles  the  capacity  region  for 
degraded  message  sets  considered  by  Korner  and  Marton.  A paper  has  been 
accepted  for  publication  In  IEEE  IT . 


1.  Capacity  Theorems  for  the  Relay  Channel  (T.  Cover,  A.  El 
Gama 1 ) 

A relay  channel  consists  of  an  Input  x^^,  relay  output 
y , channel  output  y,  and  relay  sender  x (whose  transmission  depends 
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on  the  past  symbols  y^^).  The  dependence  of  the  received  symbols  on  the 
inputs  is  p(y,yj^  . The  channel  is  assumed  to  be  memoryless. 

Hie  following  capacity  theorems  have  been  proved ; 


if  y is  a degraded  form  of  y^^ , then 


C = max  min 
pCx^.Xg) 


[i(x^,X2;Y),i(x^:Y^ 


if  y^  is  a degraded  form  of  y,  then 

C = max  max  I(X. ;Y|x  ) 

P(X^)  Xg 

If  pCy.Yj^  Ixj^fXg)  is  an  arbitrary  relay  channel  with  feedback  from  (y, 
y^^)  to  both  x^  and  x^,  then 


C 


max  min 
pCXj^.x^) 


[i(X^,X2;Y),I(X^;Y,Y^ 


j A superposition-block  Markov  encoding  shows  the  achieve- 

bility  of  C,  and  converses  are  established.  The  capacities  of  the 

■ gaussian  relay  channel  and  certain  discrete  relay  channels  have  been 

' evaluated , and  a lower  bound  on  the  capacity  of  the  general  relay  chan- 

nel  has  been  established.  A paper  has  been  submitted  to  IEEE  IT. 

t 

B.  Project  6240.  OPTIMAL  CODES  WITH  VARIABLE  WEIGHT  SYMBOLS 

\ Principal  Investigator;  J.  Gill 

I Staff:  N.  Cot 

I 1 . Objective 

1 

■ The  objectives  of  this  project  are  to  determine  procedures  for 

i producing  optimal  uniquely  decipherable  codes,  using  sjmibols  of  unequal 

costs,  and  to  study  the  properties  of  optimal  codes. 
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2 . Current  Status  of  Work 

We  have  demonstrated  that  variable-length  coding  schemes  can 
be  applied  to  several  fields  other  than  communications,  including 

• game  theory  where  they  are  associated  with  the  determi- 
nation of  optimal  gambling  strategies 

• data  allocation  in  storage  and  data  retrieval  [1]  where 
their  application  is  based  on  the  relationship  between 
double-chained  trees  and  variable-length  codes;  a file 
can  be  organized  as  a double-chained  tree  in  which  keys 
correspond  to  variable-length  codewords 

An  optimal  tree  (or  code)  corresponds  to  a minimal  average  search  time. 
If  0 < b < b . . . < b are  the  costs  of  the  symbols  of  a code  alpha- 
bet  and  A is  determined  by 


then  the  converse  to  the  source-coding  theorem  establishes  a lower  bound 
on  the  average  cost  c of  a code  (wj^,...,w^)  with  codeword  probabili- 
ties (p.,  , . . • ,P„  }, 

1 n 


where  H(p^,...,p^)  is  the  entropy  (in  bits)  of  the  probability  distri- 
bution p,  , . . . ,p  . 

1 n 

We  have  also  determined  an  upper  bound  on  the  cost  of  optimal 

codes , 


c < 


H(p^ Pj^) 

logg  + b^  + 5 


where  B < b is  a quantity  calculated  from  b , ...,b  . Using  this  up- 

"X  XX 

per  bound,  nearly  optimal  codes  can  be  realized  for  unequal  probabili- 
ties . 
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For  equiprobable  codewords , a procedure  for  producing  optimal 
codes  has  been  analyzed  [2]  and  shown  to  require  O(log^  n)  steps  to 
produce  optimal  codes  with  n codewords,  and  the  cost  of  these  codewords 
was  determined  In  Reference  3.  This  procedure  Is  being  extended  to 
nearly  equiprobable  codewords . 


References 

1.  L.  Stanfel,  "Optimal  Trees  for  a Class  of  Information  Retrieval  Prob- 
lems," Info . Stor . Retr . , 9,  1973,  pp.  43-59. 

2.  N.  Cot,  "The  Complexity  of  the  Variable-Length  Encoding  Problem," 
Sixth  Southeastern  Conference  on  Combinatorics,  Graph  Theory,  and 
Computing,  1975. 

3.  N.  Cot,  "cost  of  Optimal  Prefix  Codes,"  Johns  Hopkins  Conf.  on  System 
Sciences , Apr  1977 , 


C.  Project  6302.  PROBLEMS  IN  CONTROL 

Principal  Investigator;  G.  F.  Franklin 
Staff:  S.  Shah 


1 . Objective 


This  study  Is  concerned  with  design  problems  In  automatic  con- 
trol systems,  especially  under  conditions  of  model  uncertainty. 


2 . Current  Status  of  Work 

The  characteristics  of  the  self-tuning  regulator  have  been 
studied,  plus  related  schemes  for  simultaneous  identification  and  con- 
trol. The  goal  Is  to  develop  an  algorithm  for  adaptive  control  \dilch 
will  apply  to  multivariable  processes  that  may  be  nonmlnlmum  phase. 

The  following  guidelines  have  been  established. 


• The  structure  of  the  model  must  be  selected  to  produce 
simple  Identification  and  control  algorithms  simultane- 
ously. In  the  self -tuning  regulator,  an  ARMA  model  with 
least-square  Identification  and  minimum-response  squared 
error  constitutes  such  a structure;  however.  In  the  mul- 
tivariable case,  the  solution  is  not  yet  known. 
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• Tile  requirements  for  the  convergence  of  the  identifica- 
tion algorithm  can  be  weak.  Because  the  algorithm  is 
to  be  used  for  control , the  self -correcting  nature  of 
the  control  action  reduces  the  restrictions  on  parame- 
ter accuracy. 

• Inclusion  of  forced  delay  in  the  model  may  be  able  to 
account  for  the  nonmlnlmum-phase  characteristics  of  the 
plant . 

. Control  algorithms  must  be  developed,  idiich  are  easily 
implemented  with  varying  plant  parameters  but  still  al- 
low for  constraints  on  the  control  effort. 


D.  Project  6502.  DATA -COMPRESS ION  TECHNIQUES  AND  APPLICATIONS 

Principal  Investigator:  R.  M.  Gray 
Staff:  Y.  Linde,  Y.  Matsuyama 

1 . Objective 

The  objective  of  this  work  is  to  develop  nearly  optimal  moder- 
ately complex  data -compress ion  algorithms.  The  primary  focus  is  on  tree 
and  trellis  coding  techniques.  Common  random  process  models  and  real 
speech  data  are  considered. 

2.  Current  Status  of  Work 

this  project  is  concerned  with  the  development  and  simulation 
of  time-invariant  nonlinear  digital  filters  to  be  used  as  decoders  in  a 
tree  or  trellis  data-compresslon  system  having  a matched  tree  search 
algorithm  (such  as  a Vlterbl  or  M algorithm)  as  encoder.  A technique 
called  the  "fake  process"  approach  was  developed  wherein  one  first  at- 
tempts to  simulate  or  "fake"  a given  source  by  time-invariant  nonlinear 
digital  filtering  a memoryless  binary  source.  If  a filter  provides  a 
good  fake  process,  it  will  also  perform  well  as  the  decoder  of  a tree 
data-compresslon  system.  The  following  technique  was  developed  to  pro- 
duce the  desired  fake. 

The  register  containing  the  binary  digits  was  considered  as  a 
digital  representation  of  a real  number  between  0 and  1.  This  number 
was  passed  through  a "scrambler"  or  mapping  from  (0,1)  to  (0,1)  and  then 
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through  an  Inverse-distribution  function  mapping  corresponding  to  the 
marginal  distribution  of  the  original  source.  With  the  correct  scram- 
bler and  by  following  the  Inverse  distribution  by  a linear  digital  fil- 
ter, the  marginal  distribution  of  the  source  can  be  approximated  and  the 
power  spectral  density  matched  exactly.  This  technique  led  to  one-blt- 
per-symbol  data-compresslon  systems  that  outperformed  traditional  systems 
by  1 to  2 dB  for  gauss Ian  sources  with  and  without  memory. 

These  techniques  have  been  extended  to  "pseudo-speech"  models 
consisting  of  autoregressive  filters  driven  by  white  laplacian  noise  and 
then  to  real  speech  data  obtained  from  standard  reference  files . Ihe 
pseudo-speech  experiment  was  completed,  and  performance  of  0.5  to  1 dB 
better  than  predictive  quantization  was  obtained  (predictive  quantization 
operates  well  on  autoregressive  models) . The  system  for  compressing  real 
speech  has  been  completed , and  we  are  testing  several  adaptive  and  nonad- 
aptive  fake-process  tree  coders  on  real  speech.  Initial  results  indicate 
that  a 9 dB  SNR  is  obtainable  at  8000  blts/sec;  this  performance  Is  supe- 
rior to  most  waveform  coders  but  Inferior  to  LPC  and  APC  techniques.  Be- 
cause this  system  Is  a first  attempt,  several  possibilities  for  improve- 
ment are  being  considered. 

This  project  Is  also  supported  by  the  Air  Force  Office  of  Sci- 
entific Research. 


E.  Project  6601.  INFORMATION  THEORY  AND  INFORMATION  PROCESSING 

Principal  Investigator:  M.  E.  Heilman 
Staff:  R.  Kahn,  D.  Andelman 


1 . Objective 

The  objective  of  this  project  Is  to  Investigate  the  relation- 
ships between  Information  theory,  Information  processing,  and  complexity 
theory . 


2.  Current  Status  of  Work 

A paper  entitled  The  Wiretap  Channel  with  Feedback  Is  being 
prepared.  The  central  problem  Is  confidential  communication  over 
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multiuser  information  channels.  The  secrecy  rate  is  defined  as  the  aver- 
age equivocation  of  information  at  the  eavesdropping  terminal  (on  a per 
transmission  basis),  which  is  common  to  both  of  the  legitimate  users — in- 
cluding the  message  information  to  be  transmitted  between  them  and  any 
residual  randomness  that  may  occur  in  the  channel.  A procedure  has  been 
developed  wherein  it  is  possible  to  extract  from  the  common  Information 
a totally  secret  message  (or  key)  which  is  Independent  of  the  main  mes- 
sage. This  secret  key  can  be  exclusive-Ored  with  the  main  message,  which 
otherwise  would  be  Inadequately  protected  by  standard  wiretap  encoding. 

In  addition  to  using  this  technique  to  increase  message  equivo- 
cation when  transmitting  on  a block-by-block  basis,  it  is  also  a very  ef- 
fective procedure  when  time-sharing  between  two  modes  of  operation — one 
with  a high  key-generation  rate  and  the  other  with  a high  message-trans- 
mission rate.  A very  simple  relationship  has  been  observed  between  the 
set  of  achievable  rate-secrecy  pairs  (R,S)  and  the  corresponding  set  of 

achievable  rate-equivocation  pairs  (R,R  ) ; the  latter  is  an  alternate  but 

© 

equivalent  description  of  the  (R,d)  region.  Based  on  this  formalism,  the 
paper  will  describe  an  achievable  key-generation  rate  (assuming  public 
feedback)  and  improved  regions  for  the  independent  binary  erasure  and  in- 
dependent binary  symmetric  wiretap  channels. 

A long-held  belief  that  Slepian-Wolf  feedback  is  optimal  has 
beer  proven  false  by  counterexample.  We  are  also  extending  the  wiretap 
channel  with  feedback  to  a two-way  channel  with  eavesdropper.  In  addi- 
tion, the  feasibility  of  a convolutional  wiretap  encoder  has  been  stud- 
ied, but  no  definite  results  have  been  obtained. 

F.  Project  6701.  THE  RECURSIVE-FILTER  ADAPTIVE  LINE  ENHANCER'*’ 

Principal  Investigator:  B.  Wldrow 
Staff:  P.  Thompson 

1.  Objective 

The  adaptive  line  enhancer  (ALE) , consisting  of  a digital  filter 
adapted  to  produce  a minimum  mean-squared -error  prediction  of  a future 

^This  work  is  supported  partially  by  the  United  States  Department  of  En- 
ergy. 
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Input,  has  proven  to  be  a useful  tool  for  detecting  narrowband  signals  | 

obscured  by  wideband  noise  and  for  estimating  the  spectral  characterls-  ^ 

tics  of  signals . Earlier  research  has  provided  a good  basis  for  under-  \ 

standing  the  behavior  of  the  ALE  when  the  adaptive  filter  Is  constrained  ' 

to  an  all-zero  tapped-delay-llne  (TDL)  configuration.  It  Is  widely 
recognized  that  a recursive  filter,  containing  both  poles  and  zeros,  can 
perform  better  with  a smaller  number  of  parameters  than  an  all-zero  fil- 
ter; however,  procedures  for  adjusting  these  parameters  are  not  well  es- 
tablished. The  objective  of  this  research,  therefore.  Is  to  develop  and 
analyze  appropriate  adaptation  algorithms  for  a recursive-filter  adaptive  i 

line  enhancer. 

2.  Current  Status  of  Work 

Several  algorithms  are  being  analyzed  and  compared.  All  are 
based  on  the  well-known  relationship  between  the  optimal  filter  for  esti- 
mating a signal  In  white  noise  and  the  corresponding  Inverse  filter  that 
whitens  the  observed  process.  Ihls  relationship  makes  adaptive  filtering 
an  adaptive  Inverse-modeling  problem  encountered  In  other  such  disciplines 
as  system  Identification  and  time-series  analysis. 

The  most  general  of  the  algorithms  being  considered  Is  an  ap-  | 

proximate  maxlmum-llkellhood  scheme  similar  In  structui^  to  an  extended 
Kalman  filter.  The  gradient  of  the  mean-squared  value  of  the  Inverse- 
filter  output  (residual)  with  respect  to  the  vector  of  adjustable  filter  j 

parameters,  plus  an  Inverse  covariance  matrix.  Is  approximated  from  fil- 
tered Inputs  and  residuals  to  update  the  parameters.  The  filters  used  for 
computing  the  updates  consist  of  the  pole-only  portion  of  the  adaptive 
Inverse  filter  which  changes  during  adaptation — hence,  the  approximate 
nature  of  the  algorithm. 

The  computations  per  Iteration  required  for  this  algorithm  are 

2 

0(n  ),  where  n Is  the  number  of  poles  plus  the  number  of  zeros.  This 
can  be  reduced  to  0(n)  by  replacing  the  covariance  matrix  with  Its  trace ; 
however,  the  price  paid  for  this  reduction  Is  an  Increase  (particularly 
at  low  slgnal-to-nolse  ratios)  In  the  number  of  Iterations  required  for 
adaptation.  Compromise  methods  are  being  sought  to  produce  fast  adapta- 
tion with  minimum  computation  per  Iteration.  In  addition,  the  trade-offs 
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between  adaptation  time  and  signal-estimation  error  are  being  studied 
experimentally  and  anal3rtlcally. 

When  signals  are  sinusoidal  or  very  narrowband  (with  respect 
to  the  Nyqulst  rate),  the  most  promising  of  the  adaptation  techniques 
Is  a bootstrap  method  wherein  only  the  feedforward  parameters  are  ad- 
apted, and  the  feedback  parameters  are  slaved  to  the  feedforward  param- 
eters so  that  the  poles  (In  the  Inverse  filter)  are  slightly  closer  to 
the  origin  in  the  z-domaln  than  are  the  zeros,  thereby  facilitating  the 
desired  notch  in  the  Inverse  filter  and  peak  In  the  signal-estimation 
filter.  This  preselected  relationship  between  ptoles  and  zeros  constrains 
the  filter  away  from  the  optimal  solution  by  setting  the  minimum  attain- 
able bandwidth,  but  it  has  been  demonstrated  In  simulations  that,  for  a 
sinusoidal  signal,  the  constrained  filter  performs  as  well  or  better  than 
the  unconstrained  filter  when  algorithm  gains  are  adjusted  to  produce 
comparable  convergence  times.  It  has  also  been  observed  that  the  per- 
formance of  the  bootstrap  recursive  filter  is  comparable  to  a TDL  filter 
with  similar  bandwidth  capabilities.  Because  the  bandwidth  In  the  re- 
cursive filter  Is  not  coupled  to  the  number  of  parameters  as  it  is  in 
the  TDL,  the  recursive  filter  is  superior  in  terms  of  computations  per 
iteration.  Quantification  of  these  results  and  characterization  in  other 
respects  are  being  analyzed. 


G.  Project  7050.  STUDIES  IN  STATISTICAL  SYSTEM  THEORY 

Principal  Investigators:  T.  Kailath,  M.  Morf 
Staff:  J.  Dobbins,  D.  Lee,  B.  Levy,  J.  Newkirk, 
G.  Verghese,  E.  Verriest 

1 . Objective 


The  purpose  of  this  project  is  to  analyze  the  basic  problems 
in  statistical  system  theory. 


2 . Current  Status  of  Work 

a.  Generalized  State-Space  Systems  (T.  Kailath,  G.  Verghese, 
B.  L4vy) 

Research  has  focused  on  identifying  and  eliminating  dif- 
ficulties encountered  during  earlier  treatments  of  systems  in  the  form 
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of  Ex  s Ax  'fBu,  E Is  a square  but  singular  matrix.  Our  work 

has  led  to  a better  understanding  of  the  special  properties  of  such  sys- 
tem descriptions,  methods  for  analyzing  and  transforming  them,  and  their 
significance  In  practical  applications.  Generalized  state-space  descrip- 
tions of  the  above  form  are  common  when  the  governing  equations  are  writ- 
ten In  terms  of  the  "natural"  system  variables.  They  appear  w4ien  describ- 
ing the  Interconnection  of  subsystems  to  form  composite  systems  and.  In 
this  respect,  they  are  particularly  significant  when  displaying  the  struc- 
ture of  LCR  networks ; they  also  appear  In  the  analysis  of  singularly  per- 
turbed dynamic  systems . 

The  distinguishing  feature  of  such  systems  Is  that  they  can 
exhibit  an  Impulsive  response  to  arbitrary  Initial  conditions  at  t = 0- 
whlch  corresponds  to  natural  frequencies  at  Infinity.  The  reason  for 
this  Is  that  the  matrix  sE-A  (which  determines  natural  frequencies) 
loses  rank  not  only  for  a set  of  finite  values  of  s but  also  for  Infi- 
nite s . 

Such  behavior  does  not  appear  to  have  been  considered  be- 
fore, perhaps  because,  If  the  system  Is  In  operation  prior  to  t = 0, 
Initial  conditions  must  satisfy  the  system  equations  and  cannot  be  arbi- 
trary; the  resulting  constraints  guarantee  that  Impulsive  behavior  will 
not  occur.  Several  problems  appear,  however,  when  the  system  Is  formed 
only  at  t = 0 (as  a result  of  the  composition  of  subsystems  or  compo- 
nent failure).  As  with  LCR  networks,  the  essential  structure  may  be  such 
that,  ^ assembled  at  t = 0 with  arbitrary  Initial  conditions.  Its  be- 
havior is  Impulsive  or,  In  singularly  perturbed  structures,  the  limiting 
singular  problem  may  already  be  a generalized  system  with  arbitrary  Ini- 
tial conditions. 

For  such  s3rstems,  we  have  defined  a generalized  system 
order,  developed  concepts  of  controllability  and  observability  "at  In- 
finity," and  generated  transformations  to  obtain  equivalent  systems. 
Several  Interesting  results  can  now  be  obtained  In  this  framework. 

b.  Square-Root  Algorithms  for  Parallel  Processing  (M.  Morf, 

J.  Dobbins,  J.  Ne^lrk,  T.  Kallath) 

With  time-invariant  systems,  our  previously  reported  par- 
allel-processing results  yield  an  algorithm  that  calculates 


13 


1/2  1/2  1/2 

P (t) ,P  (2t),P  (4t),...  on  successive  iterations,  where  P is  the 

matrix  Rlccatl  variable.  Hils  rapid  method  of  calculating  the  steady- 
state  Rlccati  solution  is  valuable  «4ien  synthesizing  constant-parameter 
filters  or  controllers.  A computer  implementation  has  been  written  and 
debugged.  It  has  performed  well  on  a variety  of  test  problems,  includ- 
ing some  for  vdiich  competitive  algorithms  based  on  an  eigenvector  decom- 
position of  the  associated  Hamiltonian  system  have  failed  to  converge. 

Parallel -processing  results  have  been  obtained  when  the 
available  information  is  the  covariance  of  the  process  under  observation 
rather  than  a state-space  model  generating  the  process.  Realizing  these 
algorithms  in  square-root  form  appears  to  require  "square  roots"  of  in- 
definite matrices.  Factorizations  in  the  form  of  M = A E A'  can  be 
used,  where  E is  a signature  matrix,  and  transformations  T (such  that 
T E T'  = E)  have  been  developed  for  updating  the  factors.  The  numerical 
stability  of  these  transformations,  however,  remains  to  be  verified,  and 
Interpretation  of  the  factors  is  still  in  question.  The  standard  square- 
root  arrays  may  be  interpreted  as  arrays  expressing  certain  estimate  er- 
rors in  terms  of  orthonormal  bases.  Analogous  interpretations  of  the 
A E A'  factorizations  would  require  basis  random  variables  with  negative 
covariances.  A more  satisfactory  explanation  should  result  from  the  fol- 
lowing work. 


c.  More  General  Array  Methods  for  Least-Squares  Estimation 
(J.  Dobbins,  T.  Kailath) 

T 

Factorizations  in  the  form  of  M = UDU  (where  D is  di- 
agonal and  U is  unit  triangular)  have  attracted  considerable  attention 

because  the  computationally  expensive  scalar  square  roots  of  the  standard 

T 

square-root  factorization  M = SS  are  avoided.  From  a geometric  point 

of  view,  the  U array  expresses  estimate  errors  in  terms  of  a basis  set 

T 

having  covariance  D,  \diich  obtains  a more  compact  derivation  of  UDU 
processing  algorithms  than  is  available  in  the  literature . 

Factorizations  in  the  form  of  M =s  A E A ' , where  E may 
have  negative  diagonal  entries,  occur  in  some  parallel-processing  algo- 
rithms and  in  square-root  versions  of  the  Chandrasekhar  algorithm.  At 
first  glance,  the  above  geometric  interpretation  appears  to  imply  basis 
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vectors  with  negative  covariances;  however,  closer  study  reveals  that 
the  appropriate  vectors  are  Individually  of  positive  variance  but  cor- 
related in  such  a way  that  the  norms  of  the  A rows  are  correctly  com- 
puted via  the  E matrix.  Details  remain  to  be  Investigated,  but  It  Is 
likely  that  alternative  transformations  based  on  the  true  (positive-def- 
inite) covariance  of  the  basis  sets  may  be  more  numerically  stable  than 
"E  orthogonal"  transformations. 

d.  Realization  of  Time-Varying  Linear  Systems  (T.  Kallath, 

M.  Morf,  E.  Verrlest) 

The  concept  of  an  operational  transfer  function  was  in- 
troduced by  Kamen.  A time-varying  multiplier  can  be  Interpreted  as  an 
operator  acting  on  the  space  of  distributions.  We  consider  a ring  T 
generated  by  some  finite  set  of  such  operators  and  closed  under  differ- 
entiation. A set  of  linear  time-varying  differential  equations  can  thus 
be  represented  In  operator  form  as  Ay  = Bu,  where  the  elements  of  the 
matrices  A and  B belong  to  a noncommutatlve  polynomial  ring  T[p], 
and  p Is  the  differential  operator. 

Some  mathematical  prop>erties  of  T[p]  have  been  derived, 
particularly  a weak-dlvision  algorithm  (left  and  right)  that  leads  to  the 
notion  of  a "greatest  common  divisor"  (left  and  right)  which  Is  not  unique 
because,  generally,  T[p]  Is  not  a Euclidean  ring.  In  addition,  an  al- 
gorithm for  a least  common  multiple  (left  and  right)  has  been  obtained. 
The  existence  of  a l.c.m.  facilitates  the  "transformation"  between  LFD 
and  RED  although  problems  with  nonmonlcity  are  as  yet  unresolved. 

System-theoretic  applications  of  the  above  realizations 
Include  the  existence  of  an  operational  transformation  A ^B.  Ihese  can- 
onical realizations  also  have  a natural  representation  as  an  L/R-fractlon 
description  and  an  L/R-polynomlal  form.  The  well-known  Silverman  crite- 
rion for  Instantaneous  controllability  has  been  presented  in  a new  con- 
figuration wherein  control  Is  related  to  the  existence  of  a similarity 
transformation  that  brings  the  realization  Into  controller  form.  Condi- 
tions for  such  ring  realizability  are  being  investigated. 

The  definition  of  "mode"  has  been  extended,  and  the  prob- 
lem of  solving  a differential  equation  was  observed  to  be  equivalent  to 
the  determination  of  the  right  factors  of  a polynomial  T[p]. 
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e.  ladder  Forms  for  Least-Squares  Estimation  (M.  Morf, 

D.  Lee,  T.  Kalla th) 

Ladder  forms  are  probably  the  most  promising  canonical 
realizations  In  estimation  and  system  Identification  because  of  their 
low  computational  complexity,  stability  "by  inspection,"  and  their  re- 
lationship to  such  physical  properties  as  reflection  or  partial-corre- 
lation coefficients.  In  addition,  they  also  appear  in  scattering  and 
network  theory. 

We  have  developed  two  very  useful  recursive  ladder  forms 
for  the  exact  least-squares  linear  predictor  problem;  they  are  the  "pre- 
windowing" ladder  in  >Wiich  data  values  before  the  known  ones  are  assumed 
to  be  zero  and  the  "non-wlndowlng"  (or  covariance)  ladder  that  works  only 
within  the  available  data.  The  covariance  ladder  requires  a slight  in- 
crease In  computations,  but  the  undesirable  effects  caused  by  the  window- 
ing of  data  are  avoided . Both  ladder  forms  are  Implemented  by  exponen- 
tially weighted  error  criteria  to  treat  (slowly)  time-varying  cases,  and 
excellent  results  have  been  obtained  in  such  applications  as  speech  mod- 
eling and  high-resolution  spectral  estimation. 

Work  is  still  in  progress  to  develop  a recursive  ladder 
for  autoregressive  moving  average  (ARMA)  modeling  and  to  implement  mul- 
tichannel (multivariable)  ladder  forms. 
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A.  Project  6961.  DESCRIPTION  LANGUAGES  AND  DESIGN  TOOLS  FOR  GENERAL- 

PURPOSE  COMPUTER  ARCHITECTURES 

Principal  Investigators:  M.  J.  Flynn,  W.  M. 

vanCleemput 

Staff:  T.  Bennett,  K.  Stevens,  E.  Slutz,  W.  Cory, 

D.  Hill,  J.  Smith,  W.  Marti,  A.  von  Bechtol- 
shelm 

1 . Objective 

The  goals  of  this  study  are  to  develop  description  languages 
and  other  design  aids  to  describe  and  implement  computer  architectures 
and  to  establish  a basis  for  understanding  the  computer-design  process. 

2 . Current  Status  of  Work 

a.  The  SPRINT  Printed-Circuit  Design  System  (K.  Stevens,  T. 

Bennett,  W.  M.  vanC'eemput) 

The  objective  of  this  system  is  to  develop  an  Interactive 
computer-assisted  design  of  printed-circuit  boards.  In  the  SPRINT  sys- 
tem, the  critical  components  are  manually  placed  and  components  such  as 
14-  and  16-pin  dual-in-line  packages  are  automatically  placed.  The  in- 
terconnection routing  module  manually  routes  the  critical  connections 
and  automatically  routes  noncrltical  connections.  The  current  system  is 
limited  to  two  signal  layers;  however,  an  extension  to  multilayer  boards 
is  planned . 

The  input  to  SPRINT  is  the  Structural  Description  Language 
(SDL).  This  language  is  also  used  as  input  to  a logic  simulator,  a fault 
test-generation/simulation  system,  and  an  automatic  logic-diagram  genera- 
tion system. 

SPRINT  is  Implemented  in  MORTRAN  and  FORTRAN  IV  on  the  IBM 
370  at  SLAG  and  makes  use  of  the  Tektronix  4013  terminal.  The  SDL  com- 
piler is  implemented  in  SPITBOL  (a  SNOBOL  dialect).  The  output  is  a plot 
from  which  the  artwork  must  be  generated  manually;  however,  work  is  in 
progress  to  generate  artwork  automatically  by  means  of  a photoplotter. 
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b . Computer-Aided  Layout  of  large-scale  Integrated  Circuits 

(E.  Slutz,  W.  Marti,  W.  M.  vanCleemput) 

Although  several  systems  exist  for  the  automated  layout 
of  I£I  circuitry,  none  obtains  a layout  comparable  to  one  designed  man- 
ually. The  objective  of  this  project  Is  to  develop  and  Implement  a sys- 
tem that  uses  algorithmic  approaches  in  which  certain  decisions  are  the 
responsibility  of  the  designer.  It  is  expected  that  this  system  will 
reduce  design  time  considerably  at  no  expense  of  excessive  silicon  area. 
The  system  is  now  in  its  initial  implementation  phase,  and  the  data  base 
and  interactive  graphics  functions  are  scheduled  for  completion  in  1978. 

c.  Interactive  System  for  Design  Capture  (A.  von  Bechtol- 

sheim,  W.  M.  vanCleemput) 

In  the  current  design  system,  all  input  is  in  the  form  of 
SDL,  the  Structural  Design  Language.  Because  designers  often  prefer 
schematics  for  their  tasks,  a system  has  been  implemented  that  will  take 
as  its  input  u schematic  from  an  interactive  graphics  terminal  and  will 
output  an  SDL  description  for  further  processing  by  the  various  design- 
automation  programs.  The  current  version  is  written  in  MAINSAIL  for  the 
PDF  10  system,  and  a more  portable  version,  written  in  PASCAL,  is  being 
developed . 

d.  Implementation  of  a Digital  Design  Ixinguage  (W.  Cory,  W. 

M.  vanCleemput) 

A compiler  and  simulator  system  for  the  Digital  Design 
Language  (DDL)  is  being  implemented,  based  on  a compiler-compiler  scheme. 
The  system,  written  in  PASCAL,  should  be  reasonably  portable  and  is  sched- 
uled to  be  operational  by  the  second  quarter  of  1978. 

e.  Automated  Generation  of  Logic  Diagrams  (J.  A.  Smith,  W. 

M.  vanCleemput) 

A system  for  the  autc:natlc  generation  of  logic  diagrams, 
originally  developed  at  the  University  of  Waterloo,  Canada,  has  been  re- 
implemented and  made  compatible  with  the  SDL  input.  This  system  is  writ- 
ten in  MORTRAN  (a  structured  FORTRAN  preprocessor) . 
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f . Description  and  Simulation  of  Computer  Architectures 
(D.  Hill,  W.  M,  vanCleemput) 

The  problem  of  describing  and  simulating  a computer  S3rstem 


at  various  levels  of  abstraction  is  being  studied.  Currently,  no  adequate  I 

tools  exist  for  simulating  a design  at  multiple  levels  of  abstraction.  i 

i 

The  current  study  focuses  on  the  potential  use  of  SIMULA  for  multilevel 
simulation  of  digital  systems. 

B.  Project  7151.  COMPUTER  ARCHITECTURE 

Principal  Investigators:  E.  J.  McCluskey,  S.  S. 

Owlcki 

Staff:  D.  J.  Rossetti,  F.  W.  Terman 

1 . Objective 

The  objective  of  this  study  is  to  gain  new  insight  into  com- 
puter architecture  and  operating  systems.  We  will  investigate  design 
methods  and  programming-language  features  that  aid  in  producing  verifi- 
ably  correct  operating  systems. 

2 . Current  Status  of  Work 

A model  for  designing  and  describing  operating  systems  and 
other  concurrent  programs  has  been  developed . Programmlng-lang^uage  fea- 
tures for  implementing  such  systems  have  been  defined  and  described  by 
formal  proof  rules.  These  proof  rules  have  been  used  to  derive  verifi- 
cation methods  that  are  now  being  applied  to  sjrstem  programs  in  the  lit- 
erature . 

The  second  direction  of  this  project  is  focused  on  the  evalua- 
tion of  various  computer  architectures.  Two  methods  that  are  very  help- 
ful in  measuring  performance — a highly  accurate  modeling  technique  and 
monitoring  systems — are  being  studied. 

A set  of  trace  programs  has  been  developed  to  measure  micropro- 
cessor performance  and  to  produce  unique  Information  concerning  the  pro- 
gramming ease  of  various  microprocessors . The  tracing  is  highly  powerful 
(it  can  record  Interrupts  and  I/O  behavior)  and  almost  transparent  (ex- 
ternal hardware  reduces  significantly  the  amount  of  perturbation  intro- 
duced by  the  tracer) . 
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The  investigation  of  interleaved  memory  is  being  extended  to 
the  effects  of  request  ordering  and  to  the  most  efficient  ordering.  A 
new  model  has  been  developed  to  evaluate  the  performance  of  gracefully 
degradable  multiprocessor  systems,  based  on  such  measures  as  availabil- 
ity, processing  power,  and  proportion  of  time  in  the  degraded  mode.  The 
best  trade-off  between  computing  power  and  availability  or  between  hard- 
ware vs  software  fault  detection  can  also  be  determined . 

a.  Verifying  Concurrent  Programs  (S . S.  Owicki) 

Large  concurrent  programs,  such  as  operating  systems  or 
distributed  data  bases,  must  be  carefully  organized  if  time-dependent 
errors  are  to  be  avoided.  We  propose  a model  for  structuring  such  pro- 
grams. The  concurrent  system  consists  of  active  entities  (called  pro- 
cesses) and  passive  entities  (called  data  modules),  and  all  data  items 
must  be  local  to  a process  or  to  a data  module.  The  processes  are  or- 
dinary sequential  programs.  Data  in  modules  are  shared  by  all  processes, 
but  they  may  access  the  information  by  calling  procedures  defined  in  the 
modules.  In  addition  to  encapsulation,  modules  provide  protection  from 
time-dependent  errors  and  interprocess  synchronization.  Data  modules  can 
be  programmed  as  shared  classes  [1]. 

A system  that  fits  the  above  model  could  be  verified  by 
directly  applying  the  proof  rules  described  in  Ref.  1;  however,  this 
would  be  a very  lengthy  and  tedious  task.  To  make  verification  feasible, 
the  basic  proof  rules  must  be  used  to  obtain  more  powerful  techniques  for 
addressing  common  situations.  For  example,  we  have  developed  rules  for 
aggregating  a set  of  data  modules  and  processes  into  a single  data  mod- 
ule; proofs  of  the  system  can  then  deal  with  the  aggregate  as  a single 
entity.  A common  operating  system  technique  is  dynamic  allocation  of 
certain  resources  (such  as  I/O  devices)  among  a number  of  processes. 
During  the  time  that  a resource  is  allocated  to  a process,  it  can  be 
treated  safely  as  a local  object,  thereby  greatly  simplifying  the  pro- 
cess proof.  We  are  developing  verification  techniques  for  exploiting 
this  simplification  and  are  also  attempting  to  verify  a number  of  pub- 
lished system  programs  to  identify  other  patterns  for  which  specialized 
proof  techniques  would  be  helpful. 
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b.  Computer-System  Measurement  and  Analysis  (O.  J.  Rossetti) 

The  objective  of  this  study  is  to  gather  and  analyze  data 
concerning  computer-system  performance  to  improve  the  efficiency  of  com- 
puter architecture  and  complex  computer-program  systems.  Our  approach  is 
to  evaluate  analytically  and  empirically  past  and  future  computer  archi- 
tectures. An  essential  aspect  of  the  work  is  the  collection  of  valid 
trace  information  which  can  then  be  applied  to  specific  machine  organi- 
zations . 

Sufficient  data  are  being  gathered  to  characterize  the 
input/output  (I/O)  and  workload  environment  of  large  high-performance 
systems.  We  are  in  the  process  of  obtaining  complete  I/O  trace  infor- 
mation for  use  in  the  architecture  of  novel  storage  systems  involving 
new  addressing  schemes,  technologies,  and  system  structures.  Pour  real- 
time environments  are  being  measured,  spanning  a variety  of  operating 
systems,  processors,  and  I/O  configurations. 

The  study  has  the  following  two  phases. 


• Data  collection  at  the  various  sites:  four  large  systems 
are  being  measured — a scientific  research  center,  an  in- 
surance company,  a bank,  and  an  aerospace /electronics  man- 
ufacturer. Data  collection  Involves  extensive  operating- 
system  modifications  and  a period  (approximately  one  month) 
of  measurement  at  each  site.  This  results  in  the  most  com- 
plete set  of  such  data  yet  available  because  it  produces  a 
complete  trace  of  file  activity  over  an  extended  period  of 
time  at  both  macro-  and  microscopic  levels  of  detail.  Cur- 
rently, three  of  the  four  measurements  are  complete. 

• Analysis  of  the  data  and  their  use  in  workload  and  perfor- 
mance models:  at  the  end  of  the  measurement  phase,  we  will 
have  obtained  a large  amount  of  highly  accurate  data,  es- 
pecially regarding  I/O  and  storage -system  environments.  The 
purpose  is  to  determine  the  factors  that  comprise  a comput- 
ing workload , staying  as  far  as  possible  from  installation 
and  configuration  dependencies  and  as  close  as  possible  to 
the  "stimuli"  exerted  by  the  users  (and  their  programs)  on 
the  system.  Such  data  will  provide  a firm  foundation  for 
modeling  and  performance  evaluation  of  storage  systems  and 
new  I/O  architectures . 


i 
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Memory  Interleaving  (F.  W.  Terman) 


c . 

A series  of  models  of  Interleaved  memory  systems  has  been 
investigated  by  trace-driven  simulations.  The  basic  model  extends  the 
one  developed  by  Burnett  and  Coffman  [2]  to  the  architecture  on  the  IBM 
360/370  with  its  variable-length  InstJructlons  and  operands  [3] ; it  also 
includes  multibyte  transfers  per  memory  access.  Variations  of  this  model 
facilitate  the  study  of  the  effects  of  channel  Interference,  data-request 
reordering,  and  data-queue  emptying  at  interrupts.  By  a major  extension 
of  the  simulation  model,  the  effects  of  more  than  one  processor  accessing 
the  same  memory  can  be  analyzed. 

Memory  requests  for  the  simulations  are  obtained  from  two 
sets  of  instruction-by-instruction  trace  records.  The  first  set  traces 
the  problem  state  component  of  typical  programs  running  on  the  IBM  360/ 
370;  the  second  set  traces  samples  of  the  total  activity  of  the  CPU,  in- 
cluding both  supervisor  and  problem  states. 

The  theoretical  predictions  of  Burnett  and  Coffman  [2]  for 
the  increase  in  memory  bandwidth  c'vused  by  interleaving  agree  with  the 
simulation  results  obtained  for  the  fetching  of  instructions.  Tiie  use- 
fulness of  fetching  Instructions  in  blocks  is  limited  by  the  relatively 
high  frequency  of  branches  on  the  IBM  360/370  machines.  Consequently,  an 
active  channel  has  relatively  little  effect  (less  than  10  percent  degra- 
dation) on  the  fetching  of  instructions  from  an  interleaved  memory. 

For  the  transfer  of  operands  to  and  from  memory,  the  sim- 
ulation results  reveal  only  one-half  the  Increase  in  memory  bandwidth 
predicted  by  Burnett  and  Coffman  [2] , which  indicates  that  data  refer- 
ences on  the  IBM  360/370  are  not  random  as  was  assumed.  The  effect  of 
an  active  channel  is  again  small  because  of  the  number  of  idle  modules 
in  a typical  memory  cycle.  A larger  degradation  (10  to  20  percent)  oc- 
curs if  the  data  requests  are  forced  to  "catch  up"  with  the  Instruction 
requests  before  a successful  branch  is  executed.  This  is  a reasonable 
restriction  because  the  branch  address  cannot  be  calculated  with  cer- 
tainty until  the  values  of  the  general-purpose  registers  are  known.  On 
the  other  hand,  the  memory  bandwidth  can  be  significantly  improved  (15 
to  30  percent)  by  filling  the  data  requests  out  of  order. 
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For  high  degrees  of  Interleaving,  the  effective  memory 
bandwidth  can  be  Improved  considerably  by  the  use  of  processors;  a fac- 
tor of  1.4  Improvement  Is  typical  with  two  processors.  Additional  im- 
provement (30  to  40  percent)  can  be  obtained  by  filling  data  requests 
out  of  order;  however,  the  run  time  of  the  simulator  increases  by  as 
much  as  a factor  of  5 when  data  requests  are  reordered  because  of  the 
need  to  repeatedly  scan  the  queue  of  data  requests.  This  Increase  in 
the  run  time  of  the  simulator  Indicates  a corresponding  Increase  in  the 
operating  time  of  the  memory  manager  of  the  actual  system. 
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III.  INTEGRATED  CIRCUITS 


r 


A.  Project  5012.  IC  PROCESS  DESIGN  AND  COMPUTER  MULTILAYER 

SEMICONDUCTOR  DEVICE  ANALYSIS 

Principal  Investigator:  R.  W.  Dutton 
Staff:  D.  B.  Estreich 

1 . Objective 

The  objective  is  to  study  the  capabilities  to  predict  and  sim- 
ulate the  electrical  behavior  of  multilayer  semiconductor  devices  based 
on  layout  and  fabrication  process  modeling. 

2.  Current  Status  of  Work 

This  research  has  focused  on  distributed  multilayer  effects  on 

2 2 
I L and  CMOS  technologies . The  multilayer  effects  in  I L are  used  in  the 

design  objectives.  CMOS  multilayer  effects  can  cause  potential  latch-up 

hazards  and  must  be  avoided . 

Results  presented  recently  [1,2]  summarize  many  of  the  perfor- 

2 

mance  and  operational  limits  of  I L.  The  intrinsic  limits  related  to 
sheet  resistance  and  charge  storage  at  high  currents  can  be  modeled  di- 
rectly from  geometry  and  technology-dependent  effects.  The  extrinsic 

limits  at  low  currents  depend  explicitly  on  junction  capacitance  and  ef- 

2 

fective  current  gain.  A method  for  forming  a simplified  I L model  for 
the  extrinsic  region  is  being  pursued.  The  desired  attributes  involve 

timing  simulation  for  relatively  low-speed  control  logic  functions.  Tlie 

2 2 
facility  to  model  efficiently  I L peripheral  circuits  on  high-speed  T L 

or  EFL  chips  can  achieve  a significant  saving  in  simulation  time  «^lle 

2 

maintaining  the  essential  multilayer  effects  of  I L that  relate  to  the 
fabrication. 

The  CMOS  latch-up  study  has  focused  primarily  on  test  struc- 
ture considerations  for  detecting  potential  hazards.  Three  key  device 
issues  are 


symmetrical  path  devices  (high  aspect  ratios)  where  doping 
and  spacing  effects  dominate;  for  example,  p-well  to  source 
dimensions  and  doping  profiles  are  one  critical  combination 


• field-aided  lateral  pnp  devices  where  current-gain  can  be 

Itered  substantially  by  the  electric  field  and  basewidth 

• accurate  determination  of  sheet  resistances  such  as  p-well, 
guard  band,  and  pinched  p-well  diffusions 

This  last  item  is  important  in  defining  conditions  for  latch- 
up  initiation.  There  is  an  auxiliary  condition  to  the  unity  beta  product 
rule,  and  analysis  indicates  that  the  ratio  of  p-well  to  substrate  resis- 
tivity and  available  current  from  the  supply  are  factors  in  the  initia- 
tion. The  exact  nature  of  this  relationship  will  be  investigated. 

Tlie  principles  involved  in  the  CMOS  test  structure  and  the  mea- 
surements of  the  Stanford  test  chip  will  be  essential  in  understanding 
potential  latch-up  hazards  in  new  high  density  technology  [3]. 
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IV.  SOLID  STATE 


A.  Project  5111.  TRANSPORT  PROPERTIES  OF  AlxGai_xAs,  AlxGai-xSb  AND 

InjjGaj_j^s  SINGLE  CRYSTALS 

Principal  Investigator:  G.  L.  Pearson 
Staff:  K.  Jew 

1 . Objective 

This  project  is  a systematic  study  of  the  transport  of  free 
carriers  in  ternary  systems.  The  electronic  parameters  of  Interest  are 
free-carrier  mobility,  resistivity,  minority-carrier  diffusion  length, 
impurity-activation  energy,  and  trap-energy  levels.  The  distribution 
coefficients  of  selected  doping  impurities  as  a function  of  x are  sig- 
nificant growth  parameters.  The  results  obtained  should  prove  useful 
in  the  design  of  photo luminescent  devices,  injection  lasers,  solar  cells. 
Infrared  detectors,  and  bulk-effect  microwave  devices. 

2 . Current  Status  of  Work 

In  previous  work,  (p)Al^Ga^_^As :Ge  LPE  crystals  with  composi- 
tions throughout  the  0.0  < x < 1.0  range  were  grown  on  semi-insulating 
GaAs:Cr  substrates  and  characterized  by  Hall  measurements.  The  results 
obtained  from  these  measurements  were  generally  consistent  except  for 
crystals  of  high  aluminum  content  (x  > 0.6)  in  which  the  data  were 
scattered.  This  behavior  is  assumed  to  be  caused  by  the  difficulty  in 
doping  Al^Ga^  ^As  layers  with  x > 0.6  because  the  carrier  concentra- 
tion decreases  exponentially  with  increasing  x.  During  the  past  re- 
porting period,  (p)Al^Gaj^_^As  :Ge  was  further  investigated  by  repeating 
the  Hall  measurements  on  a separate  set  of  crystals  grown  with  higher 
dopant  concentrations. 

Before  each  growth  run,  the  Ga  melt  was  baked  out  at  800 ®C  for 
12  hours  under  flowing  H2 . The  p-type  epitaxial  layers  were  prepared 
from  Ga  melts  containing  Ge,  a GaAs  source,  and  Al.  The  dopant  concen- 
tration in  the  melt  was  increased  from  0.5  to  2.0  atomic  percent  Ge  to 
obtain  higher  carrier  concentrations  in  the  grown  layer.  The  crystals 
were  grown  in  the  temperature  interval  of  800®  to  785®C  and  at  a cooling 
rate  of  0.2®C/mln. 
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Hall  measurements  of  the  (p)Al^Gaj^  ^As:Ge  crystals  were  based 
on  the  Van  der  Pauw  technique.  Carrier  concentration,  mobility,  and  re- 
sistivity were  measured  in  the  temperature  Interval  77®  to  300 ®K.  The 
results  obtained  agree  with  previously  reported  data. 

The  hole  concentration  at  300®K  decreases  as  aluminum  content 
is  Increased . In  these  crystals  grown  from  melts  doped  with  2 percent 

Ge , the  hole  concentration  has  an  exponential  dependence  on  aluminum 

18  —3 

content  x as  p(x)  = 2.38  x 10  exp<-4.73x)  cm 

The  variation  of  carrier  concentration  with  temperature  was 
used  to  determine  the  thermal  activation  energy  of  the  Ge  acceptor.  The 
activation  energies  obtained  from  the  current  set  of  crystals  is  consis- 
tent with  = 175x  + 10  MeV.  For  samples  of  high  aluminum  content, 
the  earlier  values  for  E^  were  scattered  off  this  line;  however,  the 
present  measurements  used  crystals  grown  with  higher  dopant  concentra- 
tions, and  the  results  for  E.  vs  x are  more  consistent. 

A 

The  measured  hole  mobilities  at  300 ®K  are  in  agreement  with 
—5/2  2 

Pp(x)  = 27.8  (0.48  + 0.31x)  cm  N-sec . The  hole  mobility  decreases 

as  aluminum  content  x is  increased  because  of  the  variation  of  the  hole 

effective  mass  in  the  A1  Ga,  As  system. 

X 1-x 

Current  study  is  directed  toward  the  growth  of  (n)Al^Gaj^ 

Sn  which  will  be  characterized  by  Van  der  Pauw,  photo luminescence , I-V, 
C-V,  and  scanning  electron-microscope  techniques. 


B.  Project  5244.  STUDIES  OF  HIGH -TRANSITION  SUPERCONDUCTORS  SUCH  AS 

V3SI  AND  Nb3Sn  AND  DEPTH  PROFILING  OF  NbgSb  JOSEPH- 
SON  JUNCTIONS 

Principal  Investigator:  W.  E.  Spicer 
Staff:  I.  Lindau,  J.  Miller 


1 . Objective 

The  purpose  of  this  project  is  to  explore  the  electronic  struc- 
tures of  type-II  (hard)  superconducting  alloys,  particularly  NbgSn,  to 
determine  such  parameters  as  valence  bandwidth  uiiich  is  Important  to  the 
theoretical  models  of  superconductors.  Emphasis  is  directed  toward  the 
preparation  of  atomically  clean  Nb  surfaces,  evaporation  of  Sn  as  an  ad- 
layer,  and  annealing  to  form  Nb^Sn. 
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2.  Current  Status  of  Work 

We  have  successfully  prepared  a clean  Nb  surface  and  have  ad- 
mitted controlled  amounts  of  oxygen  to  the  system  to  observe  (via  UPS) 
the  initial  stages  of  oxidation  on  the  Nb  surface.  The  results  indicated 
that  the  photoemlsslon  spectra  at  21.2  was  not  Inconsistent  with  the  bulk 
density  of  states  calculated  by  Ho,  Louie,  Chellkowsky,  and  Cohen*^  and 
led  to  the  postulation  of  a model  for  the  formation  of  oxides  on  the 
clean  Nb  surface.  At  small  exposures  of  oxygen,  a peak  is  observed  at 
6.5  eV  from  which  is  related  to  a chemisorption  stage.  At  higher 

exposures,  a second  peak  appears  at  7.3  eV  from  which  is  associated 

with  a surface  ''bore  x = z is  a likely  candidate.  At  exposures 

up  to  saturation,  a new  peak  grows  at  8.2  eV  along  with  the  7.3  eV  peak. 
Because  the  resulting  spectra  is  similar  to  the  alr-oxldlzed  Nb  idilch  is 
known  to  consist  mainly  of  Nb20g,  the  final  phase  is  related  to  the  for- 
mation of  Nb^Og  at  the  surface.  This  research  was  presented  in  a paper 
at  the  Physics  of  Transition  Metals  Conference.* 

Recent  work  has  focused  on  the  preparation  and  construction  of 
equipment  for  new  experiments.  A more  grazing-incidence  monochrometer 
is  being  designed  and  constructed  to  improve  the  utility  and  resolution 
of  the  He  resonance  lamp.  When  this  device  is  completed,  shallow  core 
levels  of  many  materials  and  adsorbates  will  be  accessible  so  that  chem- 
ical shifts  can  be  measured.  From  these  chemical-shift  data  and  enhanced 
hlgh-resolutlon  valence-band  data,  significant  information  can  be  ob- 
tained concerning  the  chemical  state  of  the  surface.  They  can  also  be 
used  to  finger  print  different  gases  chemisorbed  on  the  metal  surface. 

The  equipment  should  be  completed  during  the  next  reporting  period,  and 
work  should  begin  on  the  formation  and  study  of  Nb^Sn  and  on  the  use  of 

sputter  Auger  techniques  to  study  Josephson  Junctions. 

I 

I 


^Phys . Rev . , B-15,  1977,  p.  1755. 

’j.  N.  Miller,  I.  Lindau,  and  W.  E.  Spicer,  Proc.  Physics  of  Trans.  Met- 
als , Toronto,  1977,  to  be  published. 
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V.  RADIOSCIENCE 


A.  Project  3606.  REAL  TIME  DIGITAL  SIGNAL  PROCESSING:  ALGORITHMS  AND 

SYSTEM  ARCHITECTURE 

Principal  Investigator:  A.  M.  Peterson 

Staff;  S.Narayan,  M.  J.  Narasimha  , A . I . Abu-el-haija 

1 . Objective 

The  objective  of  this  research  is  the  development  of  algorithms 
and  digital  system  designs  to  efficiently  and  economically  Implement 
real-time  digital  signal -processing  hardware.  The  proposed  research  is 
based  on  the  belief  that  many  signal-processing  functions  in  communication, 
radar,  sonar,  and  other  systems,  now  Implemented  by  analog  circuits,  can 
be  more  reliably  and  economically  implemented  by  digital  hardware  using 
microprocessors  and  existing  MSI/LSI  or  new  device  designs  resulting  from 
the  research . 

2 . Current  Status  of  Work 

During  the  past  two  years,  the  feasibility  and  utility  of  very 
large  digital-array  signal -processing  systems  have  been  investigated. 
This  work  was  initiated  by  a request  from  NASA  Ames  Research  Center  for 
a million-channel  spectrum  analyzer;  it  was  then  extended  to  other  ap- 
plications, such  as  processors  for  synthetic  aperture  radar  (SAR)  azimuth 
beam  forming  and  for  processing  wideband  signals  associated  with  plane- 
tary blstatlc  radar  and  radio-wave  occultation  experiments.  This  same 
processor  array  can  also  do  two-dimensional  processing  of  radar  images 
required  as  a part  of  an  ONR-supported  oceanographic  research  project 
aimed  at  sensing  ocean-surface  characteristics  by  synthetic  aperture 
radar.  These  images  are  now  obtained  from  aircraft  but  will  soon  be 
obtained  by  the  SEASAT-A  spacecraft.  This  processor  array  would  be  able 
to  perform  1000-1000  element  two-dimensional  spectral  transforms  in  a 
fraction  of  a second.  As  a near-term  application,  our  research  could 
provide  a real-time  digital  image-forming  system  for  the  SEASAT-A  syn- 
thetic radar. 
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In  the  existing  system  design  for  a million-channel  digital  ^ 

spectrum  analyzer,  digital  bandpass  filtering  transforms  a wideband  in- 
put signal  from  a single  input  channel  into  many  narrower  bandwidth  out- 

20 

put  channels.  For  example,  a single  2 Hz  a l mHz  channel  may  be  split 

20  7 13 

into  128  channels,  each  with  a (2  /2  =2  ) = 8192  Hz  bandwidth.  The 

outputs  of  this  channel-splitting  digital  bandpass-filter  bank  can  be 

7 

processed  in  2 = 128  microprocessor-based  FFT  computers,  each  producing 

13 

2 = 8192  output  channels  with  a 1 Hz  bandwidth.  The  bandpass  filters 

that  precede  the  FFT  processors  reject  interference  (>60  dB)  in  all 
channels  except  the  one  occupied  by  the  interfering  signal.  After  pas- 
sing through  the  128-channel  filter  band,  each  channel  is  Independent  of 
the  others  and  requires  identical  FFT  processing  to  obtain  the  final 
bandwidth.  The  band pass -filter  bank  that  precedes  the  parallel -process or 

array  can  be  constructed  in  a modular  form.  One  structure  that  has  been 

4 

studied  would  have  a single  input  processor  with  2 =16  output  bands  ; 

each  of  the  16  output  channels  would  be  followed  by  processors  which,  in 
turn,  provide  eight  output  channels.  The  complete  system  would  include 
128  microprocessors,  each  computing  an  8k  transform  by  performing  iden- 
tical operations  at  every  clock  cycle  but  on  different  input  data  corre- 
sponding to  the  outputs  of  the  128  bandpass  filters.  This  system  design 
is  capable  of  operating  over  a wide  range  of  input  bandwldths  up  to  e»6 
MHz. 

Recent  work  has  demonstrated  that  system  performance  can  be 
enhanced  by  using  DFT  algorithms  [1],  and  this  work  has  been  extended  at 
Stanford  [2]  and  by  Kolba  and  Parks  [3]  through  "prime  factor"  algorithms 
that  are  much  faster  than  conventional  FFT  Implementations  based  on  pow- 
ers of  two.  This  approach  is  especially  applicable  to  the  bandpass  fil- 
ters required  in  the  million-channel  analyzer.  In  particular,  a prime- 
factor  implementation  of  a 120-point  DFT  based  on  8-,  5-,  and  3-point 
DFTs  requires  only  one-third  as  many  multiples  as  a 128-point  base-2  FFT. 

Prime-factor  algorithms  should  also  facilitate  the  convolution  operations 
necessary  in  SAR  image  formation  discussed  below. 

Because  this  digital  spectrum  analyzer  is  modular,  only  three 
of  the  modules  (printed  circuit  boards)  need  be  built  to  test  the  perfor- 
mance of  the  final  system.  These  modules  are  one  of  the  128-output  FFT 
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processors , one  of  the  16  base-eight  FET  processors , and  one  base-16  pro- 
cessor (or  base  3 X 5 = 15  if  the  prime-factor  algorithm  is  used).  These 
three  building  blocks  provide  the  basis  for  a wideband  "scanning"  digital 
spectrum  analyzer  with  Impressive  capabilities.  When  the  remaining  mod- 
ules are  duplicated,  the  full  parallel  analyzer  or  any  portion  of  it  can 
be  Implemented  rapidly. 

Azimuth  processing  of  the  SEASAT-A  SAR  data  requires  beam-form- 
ing computation  and  various  doppler  and  range-walk  corrections  on  the 
4000-range  intervals  that  comprise  a 100  km  range  swath.  The  above  pro- 
cessor array  (128  processors)  would  do  the  necessary  convolution  proces- 
sing, and  a single  unit  could  handle  an  800  m range  Interval.  All  neces- 
sary azimuth  processing  on  32  range  cells  can  be  performed  in  a single  I 

processor  in  real  time.  As  a result,  the  full  100  km  range  swath  can  be 

S 

processed  in  real  time  by  120  processors.  Processing  of  approximately  1 

one-half  second  of  the  1600  pulse/sec  data  is  required  to  produce  the  1 

desired  25  m azimuth  resolution.  j 

) 

The  range  compression  of  the  SEASAT  SAR  chirped  signal  can  be  ■ 

done  digitally  in  a special  processor  that  precedes  the  multiprocessor  j 

beam-forming  system.  This  new  preprocessor  would  make  use  of  the  fast  ■ 

convolution  techniques  [3]  or  the  special  chirp-filter  design  developed 
by  Narasimha  [4]. 

The  modular  design  of  the  proposed  system  plus  I£I  results  in 
Impressive  performance  at  modest  cost.  Hie  proposed  design  is  based  on 
the  2900  family  of  bipolar  microprocessors  (4-bit  slice  architecture). 

The  recently  available  parallel  multipliers  (Monolithic  Memories  6-7558 

and  TRW  MPY  16)  and  serial  parallel  multipliers  (AMD  25L514)  also  improve 

performance . The  memory  for  data  storage  in  the  128-output  processing 

is  the  principal  cost  of  the  system.  i 
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B.  Project  4214.  INVESTIGATION  OF  UNDERSEA  COMMUNICATION  WITH  ULF/ELF 

ELECTROMAGNETIC  WAVES 

Principal  Investigator:  O.  G.  Vlllard,  Jr. 

Staff:  A.  C.  Fraser-Smith , D.  M.  Bubenik,  A.  Ho 


1 . Objective 


The  objective  of  this  project  Is  to  investigate  the  propagation 
of  ULF/ELF  electromagnetic  waves  in  the  sea  (frequencies  less  than  300  Hz) 
from  undersea  harmonic  dipole  sources,  considering  the  effects  caused  by 
the  presence  of  the  ocean  floor.  This  is  an  extension  of  earlier  work 
[1,2,3]  concerning  only  an  ocean  of  effectively  infinite  depth. 

2.  Current  Status  of  Work 


This  project  consists  of  (1)  derivations  of  the  mathematical 
expressions  for  the  fields  of  harmonic  dipole  sources  located  in  an  ocean 
of  finite  depth,  (2)  evaluation  of  these  expressions,  and  (3)  analysis  of 
the  results.  These  stages  have  been  completed,  and  a paper  entitled 
"ULF/ELF  Electromagnetic  Fields  Generated  in  a Sea  of  Finite  Depth  by  a 
Submerged  Vertically  Directed  Harmonic  Magnetic  Dipole"  by  D.  M.  Bubenik 
and  A.  C.  Fraser-Smith  is  ready  for  publication. 


ABSTRACT 

In  this  paper,  we  compute  the  undersea  quasi-static  elec- 
tromagnetic fields  of  a submerged  vertically  directed  harmonic 
magnetic  dipole  (VMD) , taking  into  account  the  presence  of  a 
sea  floor,  for  evaluating  the  relative  importance  of  the  modes 
of  propagation  of  the  fields  at  frequencies  in  the  ULF/lower 
ELF  bands  (frequencies  less  than  300  Hz).  The  source  and  the 
receiver  are  located  equidistant  from  the  surface  and  floor  of 
a sea  one  skin  depth  (6) , and  the  horizontal  separation  (r) 


between  them  is  varied  in  the  range  0.1  6 and  100  6.  When  the 
horizontal  separation  between  the  receiver  and  the  VMD  source 
is  less  than  approximately  eight  skin  depths,  each  field  com- 
ponent behaves  as  though  the  source  were  in  a sea  of  wholly 
infinite  extent.  At  larger  separations,  the  ratio  of  the  ac- 
tual field  to  the  "infinite  extent"  (or  direct)  field  ampli- 
tudes increases  rapidly  and  reaches  a value  exceeding  1036  at 
a separation  of  100  skin  depths,  indicating  that  the  fields  at 
such  distances  have  not  propagated  directly  through  the  sea 
water.  Analysis  of  the  propagation  modes  reveals  that  the  chief 
contribution  at  large  separations  (r/8  > 20)  is  the  surface 
mode,  and  the  field  amplitudes  in  this  asymptotic  regime  are 
Independent  of  sea  floor  conductivity.  At  closer  ranges  (r/6< 
10)  where  the  direct  contribution  dominates,  the  surface  and 
sea-floor  modes  are  of  equal  importance.  Between  these  limit- 
ing cases  is  a transitional  regime  where  the  first-order  sea- 
floor contribution  is  a significant,  if  not  dominant , mode . The 
changeover  between  the  transitional  and  the  asymptotic  regimes 
moves  outward  from  the  source  as  the  sea  floor  conductivity 
decreases,  and  the  field  components  in  the  asymptotic  regime 
are  characteristic  of  a quadrupole  or  higher  order  source. 
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C.  Project  4504.  TROPOSPHERIC  RADIO  PROPAGATION 

Principal  Investigator:  A.  T.  Waterman 

Staff:  R.  D.  Fleming,  T.  Torvaldson,  S.  DeLateur 


1 . Objective 


The  objective  of  this  research  is  to  investigate  atmospheric 
parameters  and  phenomena  in  the  lower  troposphere  based  on  radio  propa- 
gation in  the  S-band.  The  goals  are  to 


• measure  wind  velocity  in  the  common  volume  of  a transhorizon 
troposcatter  propagation  path  by  means  of  doppler  techniques 
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. probe  the  structure  of  winds  and  turbulence  in  the  lower 
trop>osphere  under  various  conditions  of  atmospheric  sta- 
bility and  instability 

The  plan  of  operation  involves  the  simultaneous  reception  of  3 GHz  CW 
signals  from  a transmitter  located  100  miles  away  (beyond  line-of-sight) . 
Each  of  two  narrow  beams  are  aimed  on  either  side  of  the  great-circle 
bearing  toward  the  transmitter  by  small  amounts  (1°  or  less) . If  a wind 
component  is  transverse  to  the  path  in  the  region  of  the  common  volume , 
it  is  anticipated  that  the  signal  received  on  one  beam  will  be  doppler 
shifted  to  higher  frequencies  and  that  the  signal  on  the  other  will  be 
shifted  to  lower  frequencies.  The  magnitude  of  the  doppler  difference 
should  be  proportional  to  transverse-wind  velocity.  The  two  beams  will 
be  formed  during  data  reduction  because  the  receiving  antenna  is  a sam- 
pling array  that  measures  and  stores  on  digital  tape  the  amplitudes  and 
phases  of  the  signal  received  on  12  individual  antennas.  To  obtain  high 
azimuthal  resolution,  the  antennas  are  arranged  in  a horizontal  linear 
array. 

2 . Current  Status  of  Work 

Preliminary  to  taking  data  on  the  transhorizon  path,  the  sys- 
tem was  checked  out  on  a line-of-sight  path.  The  transmitter  was  oper- 
ated from  six  locations  15  km  from  the  receiving  array.  These  six  loca- 
tions Include  two  in  line  with  the  great-circle  bearing  in  the  direction 
of  the  eventual  transhorizon  location,  two  1°  to  the  right,  and  two  1® 
to  the  left.  For  each  of  the  three  pairs,  the  separation  was  five  min- 
utes of  arc.  This  test  serves  to  verify  the  accuracy  of  the  angle-of- 
arrival  determination,  the  stability  of  the  phase  measurements,  the 
overall  operation  of  the  system,  and  the  procedures  to  be  followed  in 
the  transhorizon  experiment . 

In  preparation  for  the  data-analysls  procedures,  some  explora- 
tory computational  work  has  been  done  on  past  data  (which  used  vertical 
angle-of-arrlval  rather  than  azimuthal)  and  also  on  the  data  of  these 
line-of-sight  tests.  This  work  has  included  some  experimentation  with 
the  maximum-entropy  method  of  analysis  as  applied  spatially  to  antenna 
arrays . 
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VI.  PLASMA  PHYSICS  AND  QUANTUM  ELECTRONICS 


A.  Project  1337.  GENERATION  OF  INTENSE  MICROWAVE  RADIATION 

Principal  Investigator:  F.  W.  Crawford 
Staff:  S.  E.  Rosenthal 

1 . Objective 

The  aim  of  this  project  is  to  design  and  construct  a high-power 
magnetron  capable  of  producing  ~1  GW  ~1  ps  (~1  kJ  energy)  pulses  at  »!3 
GHz . 

2 . Current  Status  of  Work 

Work  is  continuing  on  scaling  up  in  power  the  high-power  long- 
anode  magnetron  tube  of  Boot  et  al  [1]  to  the  relativistic  electron  en- 
ergy range.  The  anode  block  has  been  designed,  a prototype  constructed, 
and  the  normal  modes  identified  by  standard  cold  testing  techniques  [2] . 
Work  is  currently  directed  towards  design  of  the  coupling-out  structure 
and  the  indirectly-heated  cathode. 

The  order  for  the  specially-constructed  high-voltage  pulse  gen- 
erator and  associated  power  supply  has  been  placed , after  detailed  con- 
siderations of  bids  by  three  different  manufacturers.  Its  cost  has  been 
covered  primarily  by  a grant  from  the  Physics  Division  of  the  NSF  and  a 
matching  grant  from  the  Stanford  University  School  of  Engineering.  It  is 
anticipated  that  the  supply  will  be  delivered  early  in  1978. 
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B.  Glnzton  Laboratory.  TWO-PHOTON  RESONANTLY  PUMPED  IR  UP-OONVERTERS 


Principal  Investigators:  S.  E.  Harris,  J.  F.  Young 
Staff;  J.  H.  Newton 

1 . Objective 

The  goals  of  this  project  are  the  development  and  extension  of 
efficient  IR  up-conversion  techniques  in  metal  vapors,  particularly  image 
up-conversion . 

2 . Current  Status  of  Work 

The  up-conversion  of  2.94  |jm  IR  images  to  4556  A has  been  dem- 
onstrated in  Cs  vaptor.  A pwwer  efficiency  of  20  percent,  with  over  10,000 

resolvable  spots,  was  obtained  using  a pump  power  of  200  kW  focused  on 
2 

an  area  of  1.1  cm  . The  Cs  cell  had  an  active  length  of  0.3  cm  and  was 

17 

operated  at  a density  of  2.0  x 10  /cc.  Images  of  a resolution  test  tar- 
get indicated  a resolution  of  better  than  12  lines/mm;  this  is  a factor 
of  2.5  of  the  diffraction  limit  of  the  optical  system  used.  The  pumping 

laser,  Nd : lanthanum-berylate , has  a natural  two-photon  coincidence  with 
2 2 

the  6s  S to  7s  S transition,  resulting  in  a simple  practical  sjrstem  with 
a number  of  potential  applications. 

During  the  last  reporting  period,  experiments  were  performed 
to  determine  the  limiting  values  of  up-conversion  efficiency  as  functions 
of  pump-power  density  and  Cs  density  to  optimize  the  system.  The  power- 
density  experiment  was  completed  and  revealed  a higher  allowed  density 
than  was  calculated.  The  Cs  density  experiment  was  only  partially  com- 
pleted because  of  a problem  with  the  Cs  cell.  Further  calculations  were 
made  on  an  extension  of  this  device,  called  a spectral  up-converter. 

The  relationship  between  efficiency  6 and  resolvable  spots 
R and  the  operating  parameters  are 


6 « 


2 2 2 
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where 

N = Cs  density 
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L s cell  length 


A a area  of  pump  beam 
P = pump  power 

To  maximize  S • R,  A can  be  made  as  small  as  possible  by  focusing 

tightly:  however,  P/A  Is  limited  by  two-photon  absorption  of  the  pump 

and  saturation  of  the  two-photon  transition.  The  calculation  of  two- 

photon  transition  probabilities  had  indicated  that  complete  saturation 

2 

should  occur  at  a power  density  of  500  kW/cm  , but  measurements  revealed 

that  the  IR  up-conversion  efficiency  was  still  increasing  quadratlcally 

2 

at  densities  of  500  kW/cm  and  increasing  near  the  3/2's  power  at  2 MW/ 

2 2 
cm  . The  increase  was  still  linear  at  3.5  MW/cm  . Because  the  tripling 

of  the  pump  (u  = 30)  ) and  the  infrared  up-converslon  process  (u  = 
s p s 

20)  + 0)  ) have  the  same  dependence  on  two-photon  absorption,  one  can 

p IR 

measure  the  tripling  efficiency  vs  pump-power  density  and  determine  the 
up-conversion  efficiency  behavior.  The  exponent  of  the  pump  power  in 
up-conversion  will  be  one  less  than  that  for  the  tripling  process.  Hils 
method  is  experimentally  easier  to  perform.  Third  harmonic  power  vs 
pump-power  density  is  plotted  in  Fig.  1. 

2 2 

At  some  Cs  density,  self-broadening  of  the  6s  S - 7s  S two-pho- 
ton transition  surpasses  the  doppler  linewidth.  Because  efficiency  is 
inversely  proportional  to  linewidth  and  the  self -broadened  linewidth  is 
proportional  to  Cs  density  N,  efficiency  at  high  densities  will  only 
increase  linearly  with  N.  Two  other  factors  related  to  increasing  N 

limit  efficiency.  First,  absorption  of  the  sum  frequency  increases  as 

2 2 0 

the  self -broadened  linewidth  of  the  6s  S - 7p  P transition  increases  and, 
secondly,  the  Cs^  molecule  density  rises  with  Cs  density  and  causes  ab- 
sorption at  the  pump  frequency.  We  have  not  yet  been  able  to  make  these 
measurements  at  sufficiently  high  densities  to  observe  such  limits.  A 
cold  spot  developed  in  the  cell  and,  even  with  maximum  heating  of  the 
side  arm,  high  enough  Cs  densities  cannot  be  achieved  to  run  the  exper- 
iment successfully.  The  cell  is  being  modified  to  correct  the  problem. 

Further  investigation  of  the  spectral  up-converter  has  yielded 
some  interesting  results.  This  process  has  the  same  two-photon  absorp- 
tion limitations  as  the  IR  up-converter  and  can  be  pumped  with  higher 
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1.079  taiMt  Denalty,  MW/cm® 


Fig.  1.  PUMP  POWER  DENSITY  SATURATION. 

pump-power  densities  than  previously  calculated.  It  is  possible  to  angle 

phasematch  the  process  over  many  coherence  lengths  and  over  a rather 

large  bandwidth  of  the  IR  by  splitting  the  pump  beam  into  two  parts  and 

recombining  them  at  the  up-converter  at  the  appropriate  angle.  Using  a 

2 18 

power  density  of  1 MW/cm  , a Cs  density  of  10  /cc,  and  a cell  length  of 
0.3  cm,  a photon  efficiency  of  greater  than  14  percent  can  be  achieved 
over  an  IR  bandwidth  of  1400  cm  . This  device  would  have  valuable  ap- 
plications for  recording  broadband  IR  spectra,  particularly  when  fast 
time  resolution  Is  required. 
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ABSTRACT 

The  time  delay  versus  frequency  of  ground  received  whistlers  pro- 
vides valuable  information  on  electron  density  and  convection  patterns 
throughout  the  magnetosphere.  In  the  past,  the  magnetospheric  environ- 
ment in  ^Miich  whistlers  propagate  has  most  often  been  modeled  by  using  a 
dipole  magnetic  field  and  field  line  plasma  density  models  with  constant 
ion-electron  temperatures.  While  the  dipole  magnetic  model  is  a good 
first  approximation  for  many  whistler  studies  of  the  magnetospheric  equa- 
torial region  out  to  ~6-7  earth  radii  (R_) , the  region  beyond  this  dis- 

£ 

tance  requires  a magne  xc  model  that  reflects  the  significant  departures 
from  a dipole  that  exist  due  to  the  interaction  of  the  geomagnetic  field 
with  the  solar  wind  and  currents  within  the  magnetosphere . Some  current 
and  future  experiments  require  an  improved  magnetic  model  to  better  re- 
late the  equatorial  positions  of  field  lines  to  their  projections  in  the 
ionosphere.  Also,  the  effect  on  whistler  studies  of  the  significant 
ionosphere-magnetosphere  plasma  temperature  differences  observed  along 
field  lines  has  not  been  explored. 

This  research  compares  whistler  analyses  based  on  the  traditional 
models  with  analyses  based  on  1)  a hybrid  quiet-time  magnetospheric  mag- 
netic field  model  that  combines  an  internal  field  model  with  a model  of 
the  average  quiet-time  contributions  to  the  internal  field  of  the  ring, 
tall,  and  magnetopause  current  systems,  and  2)  a model  of  the  plasma  den- 
sity along  field  lines  that  Includes  significant  plasma  temperature  gra- 
dients. The  major  conclusions  of  this  comparison  are: 

1)  Results  of  whistler  studies  of  the  equatorial  magnetosphere 
based  on  the  quiet-time  hybrid  magnetic  field  model  are  not  significantly 
different  from  those  based  on  a dipole  magnetic  field  model  in  the  region 
out  to  <-6-7  R^. 

2)  Beyond  6-7  R^,  whistler  measurements  of  equatorial  path  posi- 
tions, electron  densities,  and  electric  fields  based  on  a dipole  model 
are  subject  to  significant  error.  In  particular,  a)  inferred  values  of 
geocentric  equatorial  positions  of  \dilstler  propagation  paths  can  differ 
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by  0.5  to  1 R^,  b)  calculated  equatorial  electron  densities  can  differ 
by  50'5^,  and  c)  calculations  of  magnetospherlc  electric  fields  using  a 
dipole  model  can  be  a factor  of  2 or  more  smaller  than  those  based  on 
the  hybrid  model. 

3)  The  hybrid  magnetic  model  should  be  used,  Instead  of  the  dipole 
model,  In  correlative  studies  In  which  comparisons  of  ground  or  Iono- 
spheric measurements  with  magnetospherlc  measurements  are  desired.  Use 
of  the  dipole  or  Internal  magnetic  model  can  lead  to  Inconsistencies  be- 
tween the  Ionospheric  exit  positions  of  whistler  mode  signals,  as  found 
by  direction  finding  techniques,  and  the  projection  to  the  Ionosphere  of 
their  magnetospherlc  equatorial  positions,  as  found  by  whistler  analysis. 

4)  Values  of  whistler  deduced  equatorial  quantities  are  not  signif- 
icantly affected  by  temperature  profiles  In  which  the  largest  temperature 
gradients  occur  near  the  1000  km  altitude  base  of  the  field  line. 

5)  Values  of  whistler-deduced  equatorial  quantities  are  signifi- 
cantly affected  by  temperature  profiles  with  rapid  variation  of  the  tem- 
perature near  the  equator.  Near  4 R^,  If  a constant-temperature  plasma 
model  Is  used,  overestimates  of  equatorial  electron  density  by  50^  or 
more  and  underestimates  of  equatorial  path  position  by  2.5^  (0.1  R ) or 
more  can  result  from  temperature  gradients  >0.4°K/km  near  the  equator. 
These  effects  can  be  comparable  to  or  larger  than  those  resulting  from 
hybrid-dipole  magnetic  model  differences. 

6)  Whistler-deduced  electron  concentrations  at  1000  km  altitude 
can  be  increased  2 to  4 fold,  depending  on  the  nature  of  the  temperature 
gradient  at  that  altitude,  thereby  giving  better  agreement  with  certain 
low  altitude  satellite  data  on  electron  concentration. 
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NASA  NGL  05-020-008 
ONR  Nonr  N00014-75-C-1095 
and  Nonr  N00014-76-C-0689 
TR  No.  3414-3 
August  1977 

ABSTRACT 

The  study  of  wave  particle  interactions  in  the  earth's  magnetosphere 
has  advanced  markedly  in  the  last  decade.  Understanding  of  these  inter- 
actions is  important  because  of  their  possible  impact  on  the  ionosphere 
in  general  and  VLF  and  ULF  communications  in  particular,  their  control 
of  the  radiation  belt  particles  and  the  role  played  by  the  waves  in  di- 
agnostics of  the  magnetosphere.  One  important  class  of  wave  particle 
interactions  is  the  gyroresonant  interaction  of  coherent  VLF  whistler 
mode  waves  and  energetic  particles.  The  waves  involved  can  be  natural 
whistlers  or  discrete  emissions  or  signals  Injected  into  the  magneto- 
sphere from  VLF  ground  transmitters,  such  as  the  Stanford  University 
transmitter  at  Siple,  Antarctica,  and  large  scale  power  grids.  These 
coherent  waves  interact  in  the  cyclotron-resonance  mode  with  the  ener- 
getic particles  trapped  in  the  radiation  belts.  As  a result  the  waves 
grow  or  decay  in  amplitude  and  the  particles  are  perturbed  in  pitch  an- 
gle and  energy.  The  perturbations  in  pitch  angle  are  of  special  impor- 
tance in  that  they  result  in  precipitation  of  particles  into  the  iono- 
sphere. Until  recently  this  effect  has  only  been  studied  using  linear 
theory  and  analytical  techniques.  In  the  present  work,  a computer  sim- 
ulation study  of  this  Interaction  is  made  with  special  emphasis  on  com- 
puting the  wave's  effect  on  the  particles.  With  this  approach  it  is 
possible  to  obtain  a full  nonlinear  solution  of  the  equations  of  motion 
in  an  inhomogeneous  medium.  The  nonlinear  results  are  compared  with 
those  of  the  linear  theory  and  a convenient  criterion  is  presented  for 
detennlnlng  when  a complete  nonlinear  solution  is  required.  It  is  found, 
for  example,  that  in  the  case  of  equatorial  scattering  by  a 5 kHz  CW 
pulse  near  L = 4 linear  theory  begins  to  break  down  when  the  wave  am- 
plitude exceeds  3 my.  We  represent  a full  distribution  of  energetic 
particles  by  40-50,000  test  particles,  distributed  appropriately  in  phase 
space.  By  computing  the  complete  trajectories  of  all  these  particles 
the  perturbation  of  the  full  distribution  is  estimated.  The  wave  Induced 

precipitated  flux  is  computed  and  it  is  shown  that  significant  particle 
-1  2 

fluxes  (order  of  10  ergs/cm  -sec)  can  be  precipitated  into  the 
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atmosphere  by  waves  of  moderate  Intensity  (order  of  10  my).  These  fluxes 
produce  significant  perturbations  In  the  nighttime  ionosphere.  Both  the 
incoming  fluxes  and  the  ionospheric  per -urbations  appear  to  be  measurable 
by  presently  available  instruments. 
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A FEASIBILITY  STUDY 
OF  A POCKET-SIZED  GAS 
CHROMATOGRAPHIC  AIR  ANALYZER 
by  Stephen  C.  Terry 
and  John  H.  Jerman 


ABSTRACT 

The  feasibility  of  producing  a pocket-sized  air  contaminant  monitor 
based  upon  a miniature  gas  chromatograph  and  an  integrated  microcomputer 
has  been  demonstrated  and  a practical  design  for  the  instrument  is  pre- 
sented . This  pocket-sized  instrument  is  intended  to  be  used  as  a per- 
sonal exposure  monitor  which  is  unobtmislvely  carried  by  an  industrial 
worker  throughout  the  day.  The  proposed  instrument  is  fully  self-con- 
tained in  an  8 cm  X 4 cm  X 15  cm  package  weighing  approximately  0.6  kg. 

It  is  built  around  a miniature  integrated  gas  chromatograph  (GC)  consist- 
ing of  a capillary  column,  a sample  injection  valve,  and  detector.  A 
state-of-the-art  integrated  circuit  microcomputer  is  included  in  the  in- 
strument to  handle  control  and  data  processing  functions.  The  proposed 
instrument  will  be  capable  of  sampling  the  atmosphere  automatically  once 
a minute  for  eight  hours,  measuring  the  concentrations  of  up  to  10  dif- 
ferent vapors  simultaneously  to  within  10^  accuracy,  calculating  and 
storing  the  time-weighted-average  and  peak  concentrations  for  each  of 
the  gases,  measuring  concentrations  at  the  10  ppm  level,  displaying  any 
of  the  stored  concentrations  upon  demand,  and  sounding  an  alarm  when  any 
concentration  exceeds  a predetermined  value.  The  desirable  features  of 
this  instrument — its  small  size,  automatic  operation,  and  ability  to 
monitor  multiple  gases — should  make  it  a useful  analytical  tool  and  of 
significant  use  in  the  fields  of  industrial  hygiene  and  occupational 
safety. 
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INTERACTION  OF  ELECTROMAG- 
NETIC AND  ACOUSTIC  WAVES 
IN  A STOCHASTIC  ATMOSPHERE 
by  Nirdosh  Bhatnagar 


JSEP  ONR  N00014-75-C-0601 
NASA  NGL  05-020-014 
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July  1977 


ABSTRACT 

In  the  Stanford  todio  Acoustic  Sounding  System  (RASS),  an  electro- 
magnetic signal  is  made  to  scatter  from  a moving  acoustic  pulse  train. 
Under  a Bragg-scatter  condition,  maximum  electromagnetic  scattering  oc- 
curs. The  scattered  radio  signal  contains  temperature  and  wind  informa- 
tion as  a function  of  the  acoustic  pulse  position. 

In  the  theoretical  work  on  RASS  to  date,  the  effects  of  such  atmo- 
spheric parameters  as  turbulence,  humidity,  mean  temperature,  and  mean 
wind  fields  on  the  propagating  acoustic  pulse  train  have  been  ignored. 

By  neglecting  these  parameters,  the  quantitative  analyses  have  assumed 
that  the  acoustic  wavefronts  act  as  large  perfect  spherical  reflectors. 

In  this  investigation,  RASS  performance  is  assessed  in  a real  atmosphere 
where  "coherency"  of  the  acoustic  pulse  is  degraded  as  it  propagates 
vertically  into  the  lower  atmosphere.  The  only  assumption  made  is  that 
the  electromagnetic  wave  is  not  affected  by  stochastic  perturbations  in 
the  atmosphere . 

Coherency  of  vertical  acoustic-wave  propagation  is  described  through 
a perturbation-theoretic  method  and  Fejmman's  diagrammatic  technique. 

One  of  the  most  important  attributes  of  this  analysis  is  that  it  syste- 
matically and  explicitly  accounts  for  multiple  scattering  of  acoustic 
waves  in  the  presence  of  atmospheric  fluctuations  . The  coherency  results 
are  then  used  to  evaluate  the  strength  of  the  scattered  electromagnetic 
signal  from  the  acoustic  pulse  train  \diile  taking  into  account  the  pres- 
ence of  turbulence,  mean  temperature  gradients,  and  mean  wind  fields. 

It  is  concluded  that,  for  acoustic  pulses  with  carrier  frequencies 

below  a few  kilohertz  propagating  under  typical  atmospheric  conditions, 

turbulence  has  little  effect  on  the  strength  of  the  received  radio  signal 

at  heights  up  to  a few  kilometers.  This  result  implies  that  focusing  of 

RF  energy  by  the  acoustic  wavefronts  is  primarily  a function  of  sound  in- 

-2 

tensity  which  decreases  as  x , where  x denotes  altitude. 

The  effect  of  mean  vertical  wind  and  mean  temperature  on  the  strength 
of  the  received  signal  is  also  demonstrated  to  be  insignificant.  Mean 
horizontal  winds,  however,  shift  the  focus  of  the  reflected  electromag- 
netic energy  from  its  origin,  resulting  in  a decrease  in  received  signal 
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level  when  a monostatic  RF  system  is  used.  For  a bistatic  radar  config- 
uration with  space-diversified  receiving  antennas,  the  shifting  of  the 
acoustic  pulse  makes  possible  the  remote  measurement  of  the  horizontal 
wind  component . 


SEL-77-031 


A SIIilCXJN  ABSOLUTE  PRESSURE 
TRANSDUCER  FOR  BIOMEDICAL 
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ABSTRACT 

Accurate  measurement  of  pressure  is  important  for  many  diagnostic, 
surgical  and  patient  care  procedures  and  also  for  the  control  and  inter- 
pretation of  many  types  of  physiological  research  experiments.  To  accu- 
rately measure  these  pressures  a transducer  that  was  smaller,  more  stable 
and  less  expensive  than  currently  available  units  was  desired. 

A silicon  absolute  pressure  transducer  has  been  developed  which 
meets  this  need.  The  miniature  transducer  (1.0  mm  X 1.25  mm  X 0.4  mm) 
can  be  totally  batch  fabricated  using  Integrated  circuit  processing  tech- 
niques, thus  yielding  a low  cost  device.  The  pressure  transducer  consists 
of  a silicon  chip,  containing  a silicon  diaphragm  with  p-type  resistors 
diffused  into  it,  which  is  hermetically  sealed  to  a glass  cap.  An  anodic 
bonding  process  is  used  for  this  sealing  and  results  in  a very  stable 
bond.  The  use  of  inherently  stable  materials  resulted  in  a pressure 
transducer  that  was  very  stable  with  time.  The  long-term  drift  has  been 
measured  as  being  less  than  1 torr  per  month. 

Analysis  of  the  stress  patterns  on  the  silicon  diaphragm  combined 
with  analysis  of  the  piezoresistive  effect  resulted  in  positioning  and 
orienting  the  p-type  resistors  so  as  to  maximize  the  transducer  pressure 
sensitivity.  The  sensitivity  of  this  device  (17  uV/V^/torr)  is  very 
high,  especially  considering  its  small  size. 

An  analysis  of  the  effect  of  resistor  Junction  depth  on  pressure 
sensitivity  shows  that  this  effect  can  be  substantial  for  miniature 
transducers.  The  sensitivity  decreases  as  a linear  function  of  the  ratio 
of  Junction  depth  to  diaphragm  thickness.  For  maximum  sensitivity  the 
resistor  Junction  depth  should  be  less  than  10^  of  the  diaphragm  thick- 
ness . 

The  silicon  absolute  pressure  transducer  has  been  tested  both  in 
vitro  and  in  vivo.  Results  of  extensive  in  vitro  testing  have  shown  that 
the  transducer  meets  or  exceeds  all  the  criteria  for  a pressure  trans- 
ducer used  for  chronic  implants.  Results  of  acute  in  vivo  testing  have 
demonstrated  that  the  pressure  transducer,  packaged  for  biomedical  use, 
could  survive  the  rigors  of  surgical  implantation  and  give  accurate  pres- 
sure recordings.  As  a result  of  these  experiments,  it  is  felt  that  the 
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silicon  absolute  pressure  transducer  will  find  widespread  use  In  many 
biomedical  applications. 
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K ITCH IP — A SUPERIOR 
LOW-POWER  MONOLITHIC 
DEVICE  ARRAY 
by  Stephen  Richard  Combs 


ABSTRACT 


Custom  integrated  circuits  (ICs),  compared  to  discrete  component 
designs,  substantially  reduce  size  and  power  consumption  \diile  improving 
performance,  reliability,  and  device  matching.  They  require  sizeable 
development  time  and  expense,  however,  and,  as  a result,  are  not  feasible 
for  small-scale  low-production  volume  applications.  These  shortcomings 
are  alleviated  by  a new  quick  cost-effective  monolithic  array  of  uncon- 
nected transistors,  resistors,  and  MOS  capacitors  called  a KITCHIP  that 
can  be  transformed  into  virtually  any  simple  digital  or  analog  circuit 
configuration  with  a single  customizing  aluminum  interconnect  pattern. 

Two  KITCHIP  arrays  described  in  this  report  are  tailored  specifically 
for  low-power  biomedical  telemetry  systems  for  which  no  commercially 
available  ICs  exist.  A bipolar  KITCHIP  has  been  applied  in  25  different 
circuit  designs,  resulting  in  a bidirectional  CW  blood  flowmeter,  a six- 
channel  physiological  parameter  telemetry  system,  micropower  command  re- 
ceivers, and  an  ingestlble  gastrointestinal  telemetry  pill. 

Essential  device  characteristics  for  currents  between  1 pA  and  1 mA 
typical  of  these  battery-powered  systems  are  high  transconductance,  good 
low-current  cutoff  frequency,  and  minimal  noise.  A comparison  of  bipolar 
npn,  JFET,  NMOS , DMOS , and  VMOS  transistor  performance  under  common  pho- 
tolithographic constraints  Identifies  bipolar  and  VMOS  devices  as  best 
suited  to  biotelemetry  applications.  Based  on  this  conclusion,  an  ad- 
vanced integrated-circuit  technology  featuring  two  novel  complementary 
blmodal  transistors,  whose  MOS  or  bipolar  operation  mode  is  determined 
by  the  metal  Interconnection  pattern  alone,  is  developed  and  Implemented 
in  a second  monolithic  KITCHIP  array.  Electrical  characteristics  of 
these  transistors  and  related  KITCHIP  devices  are  analyzed  in  terms  of 
the  physical  parameters  of  their  vertical  impurity  profiles,  and  an  im- 
purity distribution  optimized  for  biomedical  telemetry  applications  is 
derived  on  a layer  by  layer  basis.  The  resulting  n-type  bimodal  tran- 
sistor features  alp  double-diffused  channel  length,  a 1 V threshold 
voltage,  and  bipolar  current  gain  of  200.  P-type  transistor  character- 
istics are  only  slightly  Inferior. 
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An  eight-mask  fabrication  sequence  presented  for  realization  of 
these  profiles  employs  ion-implantation  for  excellent  control  of  electri- 
cal performance,  low-temperature  masking  oxidations  to  minimize  unwanted 
diffusion,  and  self-aligning  techniques  to  reduce  parasitic  capacitance. 
Approximate  solutions  for  processing  parameters  under  arbitrary  perfor- 
mance criteria  are  deduced  from  given  grapihical  relationships.  Novel 
bimodal  KITCHIP  device  designs  and  die  layout  which  utilize  computer- 
aided  artwork  for  accurate  geometry  matching,  incorporate  techniques  for 
enhancing  flexibility  and  versatility  such  as  scaling  of  geometry-depen- 
dent transistor  parameters  and  pre integration  of  key  circuit  configura- 
tions. With  the  bimodal  KITCHIP,  the  operation  mode  of  the  device,  its 
transconductance,  the  ratio  of  MOS  to  bipolar  parts,  and  their  relative 
placement  across  the  die  are  all  uniquely  determined  by  the  metal  inter- 
connection pattern  alone. 

Each  4 mm  square  bimodal  KITCHIP  die  contains  nearly  300  transis- 
tors, over  14.5  Mf!  of  ion-implanted  resistors,  120  pF  of  MOS  capacitance, 
and  three  device  interconnect  levels.  It  has  a unique  ability  to  inte- 
grate circuit  designs  combining  linear,  quasi-linear , digital,  and  ana- 
log-multiplex  functions  on  a single  IC  in  a quick  reliable  cost-effective 
manner . 
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ABSTRACT 

One-way  functions  are  functions  which  are  easily  computed  but  dif- 
ficult to  invert.  These  functions  have  several  uses  in  cryptography. 
They  can  be  used  to  provide  a user  authentication  system  that  is  secure 
against  theft  of  the  password  list,  and  they  are  necessary  to  the  exis- 
tence of  cryptographic  systems  which  must  withstand  a known  plaintext 
attack . 

Exponentiation  over  finite  fields  is  a candidate  one-way  function, 

and  in  addition  can  be  used  to  provide  a cr3rptosystem  vdiich  is  secure 

against  a known  plaintext  attack  if  and  only  if  the  inverse  function, 

logarithms  (indices)  over  finite  fields,  is  difficult  to  compute.  For 

the  finite  field  GF(p),  where  p is  prime,  the  best  previously  published 

1/2 

algorithm  requires  0(p  ) complexity  in  both  time  and  storage  to  com- 

pute a logarithm.  Two  improved  algorithms  for  computing  logarithms  over 
GF(p)  are  presented.  They  restrict  the  choice  of  p if  exponentiation 
over  GF(p)  is  to  be  a one-way  function  and  are  useful  in  applications 

which  require  the  computation  of  logarithms  over  GF(p) . One  of  these 

2 

algorithms  requires  only  0(log  p)  complexity  in  time  and  storage  for 
certain  primes  p. 

One  of  the  algorithms  for  computing  logarithms  over  finite  fields 
requires  the  solution  of  s3rstems  of  linear  congruences  mod  N where  N is 
a composite  number.  An  algorithm  is  developed  for  solving  these  congru- 
ences. Upper  and  lower  bounds  are  derived  for  the  fraction  of  matrices 
which  have  inverses  mod  N (1  .e . systems  of  congruences  which  have  unique 
solutions) . 

Another  candidate  one-way  function  is  the  knapsack  problem,  taken 
from  the  class  of  NP-complete  problems.  In  order  for  this  function  to 
be  difficult  to  invert,  restrictions  must  be  placed  on  the  choice  of  the 
knapsack.  Several  such  restrictions  are  derived. 

Finally,  several  methods  for  cryptanalyzing  simple  rotor  (cipher) 
machines  are  presented  for  known  plaintext  and  ciphertext-only  attacks. 
Equivalence  classes  of  rotors  are  also  investigated. 
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A SIMULATOR  FOR  THE 
EVALUATION  OF  DIGITAL 
SYSTEM  RELIABILITY 
by  Peter  Alan  Thompson 

1 

ABSTRACT  ' 


This  report  describes  a simulation  package  designed  to  evaluate  the  | 

reliability  of  digital  systems.  The  simulator  can  be  used  to  model  many 
different  types  of  systems,  at  varying  levels  of  detail.  The  user  is 
given  much  freedom  to  use  the  elements  of  the  model  in  the  way  best 

1 

suited  to  simulating  the  operation  of  a s3rstem  in  the  presence  of  faults. 

The  simulation  package  then  generates  random  faults  in  the  model  and  uses 
a Monte  Carlo  analysis  to  obtain  curves  of  reliability.  Three  examples  ! 

i 

are  given  of  simulations  of  digital  systems  which  have  redundancy.  The 

i 

difference  between  this  type  of  simulation  and  other  simulation  tech-  ^ 

nlques  is  discussed.  j 
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ABSTRACT 


This  repKjrt  presents  the  results  of  an  investigation  of  a proposed 
method  for  the  controlled  artificial  generation  of  ultra -low-frequency 
(ULF)  hydromagnetic  waves,  primarily  of  class  Pc  1 (0.2  to  5 Hz),  in  the 
ionosphere  and  magnetosphere.  The  basis  of  this  method,  which  is  called 
the  "peninsula  method"  (a  second  possible  method,  the  "VLF  method,"  is 
discussed  in  a companion  rep>ort),  is  the  passage  of  a ULF-modulated  elec- 
tric current  around  a relatively  nonconducting  peninsula  in  the  sea  or 
in  a large  saline  lake  to  form  a ULF  current  loop  that  produces  a ULF 
magnetic  field  in  the  lower  ionosphere.  Provided  the  amplitude  of  the 
ULF  magnetic  field  fluctuations  is  sufficiently  large,  i.e.,  provided 

the  maximum  magnetic  moment  of  the  peninsula  current  loop  is  greater 
13  2 

than  about  10  Am  , it  is  predicted  theoretically  that  ULF  hydromagnetic 
waves  can  be  generated  in  a disturbed  region  of  the  lower  ionosphere 
above  the  peninsula.  These  waves  can  then  propagate  away  to  large  dis- 
tances in  the  ionosphere  and  magnetosphere . 

The  peninsula  method  is  a version  of  a particular  class  of  ULF  wave 
generation  methods  based  on  the  use  of  large  ground-based  ULF  current 
systems.  Compared  with  other  possible  methods  of  generation,  these  meth- 
ods appiear  to  have  the  advantage  of  reliability  and  versatility.  How- 
ever, both  the  construction  costs  and  the  power  requirements  for  these 
systems  are  large.  The  peninsula  method  is  particularly  attractive  be- 
cause it  would  minimize  these  latter  disadvantages. 

Experiments  conducted  in  1975  and  1976  with  a small  peninsula  on 
Chappaqulddick  Island,  Massachusetts,  show  that  the  sea  (or  salt)  water 
surrounding  a peninsula  can  Indeed  function  as  a conducting  loop  and  that 
this  loop  can  be  used  to  produce  ULF  magnetic  fields  above  the  peninsula. 
A modeling  study  based  on  the  results  of  these  experiments  indicates  that 
the  peninsula  method  is  very  efficient;  the  magnetic  field  produced  at 
E region  height  by  an  electric  current  flowing  through  the  sea  water  sur- 
rounding the  peninsula  can  be  up  to  49  times  larger  than  the  magnetic 
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field  that  would  be  generated  by  the  same  current  flowing  through  a wire 
loop  laid  along  the  shoreline  of  the  peninsula.  In  addition,  because  of 
the  low  resistance  of  the  sea  vater  path,  the  power  required  to  drive 
the  current  through  the  sea  water  around  the  peninsula  can  be  more  than 
an  order  of  magnitude  smaller  than  the  power  required  to  drive  the  same 
current  through  the  wire  loop.  Thus  the  magnetic  field  produced  per 
unit  of  Input  power  is  substantially  higher  for  the  peninsula  current 
loop  than  it  Is  for  a large  horizontal  wire  loop  on  the  ground.  The 
estimated  cost  of  constructing  a peninsula  ULF  generator  Is  found  to  be 
over  an  order  of  magnitude  less  than  the  estimated  cost  of  a horizontal 
wire  loop  system  of  similar  capability.  We  conclude  that  the  peninsula 
method  of  ULF  wave  generation  is  feasible  and  that  further  experiments, 
particularly  a full-scale  ULF  generation  experiment,  are  desirable. 


SEL-77-041 


Iv 


r 


! 


¥ 


STIMULATION  OF  ULF 
GEOMAGNETIC  PUI^ATIONS 
BY  CONTROLLED  VLF 
TRANSMISSIONS  INTO 
THE  MAGNETOSPHERE 
by  A.  C.  Fraser-Smith , 
R.  A.  Helllwell,  T.  F. 
Bell,  T.  L.  Crystal, 
and  B.  Dingle 


Ahi-A  Order  No.  1733 
through  ONR  Contract 
N00014-75-C-1095 
Task  Area  NR  089-121 
and  NSF  through  OPP 
Grant  GV-28840X  and 
ARS  Grant  GA-32590X 
TR  No.  4207-8 
December  1977 


ABSTRACT 


This  report  presents  the  results  of  an  investigation  of  a proposed 
method  for  the  controlled  artificial  generation  of  ultra-low-frequency 
(ULF)  hydromagnetic  waves  of  class  Pc  1 (0.2  to  5 Hz)  in  the  ionosphere 
and  magnetosphere.  In  this  method,  which  is  called  the  "VLF  method,"  a 
large  ground-based  very-low-frequency  (VLF)  transmitter  is  used  to  stim- 
ulate the  ULF  waves  by  injecting  pulses  of  VLF  waves  into  the  magneto- 
sphere. A second  fxjssible  method  of  ULF  wave  generation,  the  "peninsula 
method,"  is  discussed  in  a companion  rep)ort . In  the  VLF  method,  the 
basic  stimulation  frequency  is  the  pulse  repetition  (or  modulation)  fre- 
quency, which  is  chosen  to  fall  in  the  ULF  range.  Stimulation  of  ULF 
waves  may  also  occur  at  any  of  the  harmonics  of  this  basic  frequency, 
dep>ending  on  the  magnetospheric  and  ionospheric  conditions  at  the  time, 
or  it  may  conceivably  occur  at  a nonharmonic  frequency  if  the  conditions 
are  particularly  favorable  and  the  generation  process  produces  a high 
level  of  broadband  ULF  hydromagnetic  noise. 

Theoretical  and  experimental  progress  with  the  VLF  method  has  been 
encouraging.  Theoretical  progress  includes  (1)  computer  simulation  of 
the  cyclotron  resonance  interaction,  which  is  believed  to  play  a key 
role  in  both  the  ULF  generation  mechanism  and  the  amplification  of  VLF 
and  ULF  waves,  (2)  completion  of  a theory  for  ULF  wave  generation  based 
on  the  repetitive  precipitation  of  energetic  electrons  into  the  iono- 
sphere by  periodic  VLF  pulses  from  a ground-based  transmitter,  and  (3) 
development  of  a method  for  computing  the  fluxes  of  energetic  electrons 
precipitated  by  VLF  waves  in  the  magnetosphere . 

Experimental  progress  Includes  the  completion  of  several  long-dura- 
tion VLF  transmission  experiments  using  the  100  kW  transmitter  at  Slple 
Station,  Antarctica.  Although  it  does  not  yet  appear  possible  to  stimu- 
late ULF  waves  on  demand  using  this  VLF  transmitter,  our  experiments 
suggest  that  VLF  transmissions  from  Siple  could  alter  the  occurrences 
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of  Pc  1 pulsations  at  Roberval,  Quebec,  which  is  geomagnetically  conju- 
gate to  Siple  Station.  Because  the  Pc  1 pulsation  events  that  occur  at 
Roberval  can  also  probably  be  observed  on  many  occasions  over  much  of 
North  America  (and,  by  inference,  over  the  equivalent  conjugate  area  in 
the  Southern  Hemisphere),  our  experiments  suggest  that  VLF  transmissions 
into  the  magnetosphere  from  a single  appropriately-located  transmitter 
can  Influence  ULF  activity  over  a large  area  of  the  earth's  surface. 

Thus,  further  experiments  appear  to  be  justified. 

Combining  the  theoretical  and  experimental  results  obtained  during 
this  research,  it  is  suggested  that  naturally-occurring  repetitive  VLF 
activity  can  stimulate  Pc  1 pulsation  events,  and  it  is  further  suggested 
that  such  VLF  activity  may  be  a major  source  of  stimulation  for  Pc  1 
pulsations.  Thus,  future  experiments  on  ULF  wave  generation  with  ground- 
based  VLF  transmitters  would  probably  benefit  greatly  if  they  were  com- 
bined with  a program  of  simultaneous  observations  of  naturally-occurring 
VLF  and  ULF  activity. 
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ABSTRACT 


! In  medical  diagnostics,  it  is  extremely  useful  to  have  available  an 

! Instrument  that  can  produce  real-time  images  of  organs  in  the  human  body. 

Such  a device  enables  detailed  investigation  of  heart  motion  and  facili- 
tates transducer  placement  for  maximum  image  quality.  Low-power  ultra- 
! sonic  waves  used  as  the  sensing  radiation  make  possible  noninvasive  ex- 

amination without  any  known  patient  risk. 

The  principal  objective  of  this  research  is  the  development  of  an 
electronically  scanned  phased-array  ultrasonic  imaging  system  based  on 
CCD  delay  elements.  Charge-coupled  devices  can  achieve  variable  analog- 
signal  delays  with  large  time -bandwidth  products  on  a small  silicon  chip, 

* and  they  are  the  key  factor  in  the  potential  realization  of  a compact 

^ • flexible  imaging  system.  The  phased  array  contains  31  piezoelectric 

transducers  and  a separate  transmitter  and  receiver  for  each  element. 

..  Bh^ery  receiver  incorporates  a CCD  to  synthesize  electronically  the  equiv- 

alent of  a lens  focused  on  the  desired  object  point.  Additional  circuitry 
sums  the  element  outputs,  controls  array  scanning,  and  generates  the  dis- 
play. The  results  are  the  ssmthesis  of  an  ultrasonic  beam  pattern  sweep- 
ing over  a 60°  planar  sector  in  the  body  and  an  oscilloscope  display  of 
the  received  ultrasonic  echoes  from  a cross  section  of  the  human  organ. 

Doppler  beam  distortion  caused  by  dynamic  focusing  and  phase-modu- 
lation sidebands  resulting  from  nonuniform  clocking  of  a CCD  have  been 
analyzed.  It  has  been  found  that  the  donpler  effect  increases  system 
resolution  slightly  at  near  ranges  and  lowers  the  average  echo  frequency. 

An  experimental  system  operating  with  25  active  channels  has  been 
tested,  and  measured  beam  patterns  are  presented.  It  utilizes  continu- 
|!  ous  variations  of  the  CCD  delay  with  clock  frequency  to  produce  dynamic 

■;  focusing.  The  multiplicative  algorithm,  optimized  for  single-input/ 

’ single-output  CCDs,  obtains  real-time  scan  generation  and  control. 
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ABSTRACT 

Integrated  circuit  fabrication  technology  has  permitted  the  devel- 
opment of  an  accelerometer  weighing  less  than  0.02  gm,  in  a 2 X3  xO.6 
mm  package.  The  accelerometer  is  intended  for  biomedical  applications 
such  as  measurement  of  heart  wall  motion  and  intrauterine  motion,  where 
extremely  small  mass  is  necessary  to  avoid  measurement  artifact.  In 
addition,  such  applications  require  a transducer  which  will  detect  ac- 
celerations down  to  0.01  g over  a 100  Hertz  bandwidth,  with  an  upper 
acceleration  limit  of  50  g. 

The  accelerometer  is  a glass-silicon-glass  sandwich  whose  active 
element  is  a very  thin  (15  pm)  cantilevered  beam  of  silicon.  A silicon 
or  gold  mass  is  mounted  at  the  free  end  of  the  beam.  A 200pm  thick  sil- 
icon supporting  rim  surrounds  the  beam  and  the  mass,  providing  the  built- 
in  end  condition  for  the  cantilever  and  space  for  contacts.  A resistor 
diffused  into  the  top  surface  of  the  beam  changes  its  value  with  accel- 
eration due  to  the  stress  induced  in  the  beam.  A second  resistor,  placed 
in  an  unstressed  region,  is  used  for  temperature  compensation.  Wells 
etched  into  the  glass  caps,  which  are  hermetically  sealed  to  the  sup- 
porting rims  using  anodic  bonding,  form  a cavity  completely  enclosing 
and  protecting  the  beam.  Fabrication  is  a batch  process  utilizing  stan- 
dard IC  photolithographic  and  diffusion  techniques,  plus  anisotropic 
etching  of  silicon,  to  shape  the  accelerometer. 

Experimental  results  agree  closely  with  the  theoretically  derived 
dep)endence  of  output  on  acceleration.  The  observed  resistance  change 
is  proportional  to  the  applied  acceleration,  and  depends  on  the  length, 
thickness,  and  width  of  the  beam,  and  on  the  magnitude  and  placement  of 
the  mass.  The  dependence  on  beam  geometry  allows  the  sensitivity  to  be 
readily  varied  over  several  orders  of  magnitude  yet  remain  tightly  con- 
trolled. Increasing  the  sensitivity  is  shown  to  decrease  proportionally 
the  maximum  acceleration  the  accelerometer  can  withstand.  The  resonant 
frequency  also  drops,  reducing  the  useful  bandwidth.  These  limitations 
are,  however,  equally  applicable  to  other  strain  gauge  accelerometer  de- 
signs. The  miniature  accelerometer's  frequency  response  is  basically  that 
of  an  ideal  tvwj-pole  system,  with  the  resonance  typically  between  500  and 
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2500  Hertz.  Although  the  damping  factor  is  only  0.005  with  air  in  the 
beam  cavity,  0.7  critical  damping  can  be  achieved  with  the  addition  of 
a fluid . 

_5 

Accelerometers  with  sensitivities  varying  from  ZsR/R  = 5 x 10  to 
-3 

2 X 10  per  g have  been  fabricated,  so  accelerations  less  than  0.001  g 
can  be  detected.  The  corresponding  operating  ranges  are  ± 200 g to  ±30 
g.  The  miniature  accelerometers  compare  very  well  with  the  small  strain 
gauge  accelerometers  available  commercially,  while  providing  more  than 
an  order  of  magnitude  reduction  in  volume  and  mass.  The  small  size  of 
this  accelerometer,  coupled  with  its  performance  and  the  low  cost  poten- 
tial of  batch  fabrication,  makes  it  extremely  attractive  for  biomedical 
use . 
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ABSTRACT 


For  the  past  several  years  we  have  carried  out  a series  of  experi- 
mental and  theoretical  investigations  of  domain  wall  motions  in  yttrium 
orthoferrite,  which  is  a low-loss,  high  mobility  canted -antiferrimagnetic 
material.  Using  modified  versions  of  the  Sixtus  and  Tonks  transit  time 
technique,  we  measured  the  domain-wall  velocity  versus  applied  H field 
behavior  of  Bloch,  Neel,  and  Head-to-Head  walls  in  bulk  single -crystal 
bars  of  YPeO„  for  values  of  applied  H field  as  high  as  100  o.e.  wall  ve- 
locities  as  high  as  10  cm/sec  and  at  temperatures  varying  from  200°K  to 
600 °K. 

Experimental  results  on  single  crystal  bars  1 mm  in  cross  section 
and  50  mm  long  show  that  the  velocity  vs  field  characteristics  for  all 
three  types  of  wall  display  a well-defined  low  field  region  where  veloc- 
ity varies  linearly  with  H,  hence  characterizable  by  a differential  mo- 

3 

billty  p.  We  find  p to  be  approximately  6 X 10  cm/sec-o.e.  for  both 

/s  4 

Bloch  and  Neel  walls  but  appears  to  be  4.2  x 10  cm/sec  o.e.  for  the 

Head-to-Head  wall,  several  times  larger  than  for  Bloch  or  Neel  walls. 
Existence  of  a linear  V-H  region  suggests  that  the  dominant  damping  pro- 
cess in  the  bulk  are  linear  ones.  When  temperature  is  increased,  mobil- 
ity, p,  for  all  3 types  of  wall  is  found  to  diminish  monotonlcally . 

In  addition,  velocity,  v,  for  all  3 types  of  wall  saturates  at  a 

5 

common  value  of  V =4.2x10  cm/sec  and  remains  pinned  to  V until 
s s 

H applied  becomes  very  large  (>750  o.e.).  This  phenomena  was  then  dis- 
covered to  be  a result  of  strong  interaction  between  the  moving  wall  and 
the  acoustic  waves  in  the  solid  when  velocity  of  the  wall  approaches  the 
propagation  speed  of  the  acoustic  modes. 

In  the  process  of  analyzing  our  experimental  results,  we  carried 
out  a systematic  computation  of  spin-wave  spectra  for  and  found 

that,  for  the  case  of  the  uniformly  magnetized  material,  two  tsrpes  of 
magnon  modes  are  present;  a rocking  mode  (a)  and  a twisting  mode  (7)* 
When  any  of  the  3 types  of  domain  walls  are  present,  each  magnon  mode 
splits  into  two  branches,  a free  spin-wave  branch  of  high  excitation 
energies  and  a bounded  wall-excitation  branch  of  relatively  lower  ener- 
gies. Similarities  of  the  magnon  spectra  in  the  presence  of  these  3 
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different  types  of  wall  lead  us  to  conclude  that  the  Gilbert  damping 
parameter  a which  describes  the  magnitude  of  energy  loss  In  wall 
motion  would  be  very  much  the  same  for  all  these  3 cases. 


Magnetostatic  energy  analysis  then  reveals  that  the  Head-to-Head 
wall  Is  actually  slanted  at  an  angle  Q = 78°  in  the  sample  bar;  and 
such  a slanting  can  be  shown  to  produce  exceptionally  high  apparent  ve- 
locities and  mobilities.  We  compensated  for  such  amplification  effects 
from  slanting  and  found  that  the  Head-to-Head  wall  Is  actually  very  sim- 
ilar to  Bloch  and  Neel  Walls  In  Its  low-field  behavior.  Magnetostatic 
analysis  also  shows  that,  for  the  particular  size  and  shape  of  our  sam- 
ple, the  Bloch  and  Neel  walls  would  experience  an  effective  Internal 
field  of  magnetostatic  origin,  propelling  them  forward  even  when  the 
external  field  Is  zero.  This  leads  to  an  anomalous  zero-field  propaga- 
tion of  Bloch  and  Neel  walls  which  was  also  observed  experimentally.  We 
suppressed  these  magnetostatic  effects  by  changing  the  shape  and  size  of 
our  sample  and  obtained  decisive  V-H  data  for  the  Neel  wall  from  a tem- 
perature of  145°K  to  390°K.  These  results,  from  earlier  analysis,  are 
expected  to  be  very  close  to  those  for  Bloch  and  Head-to-Head  wall  also. 

In  particular,  wall  mobility  at  room  temperature  Is 

4 

H s 1.2  X 10  cm/sec-o.e. 

Corresponding  to  a Bloch  relaxation  time  of 

T = 1.8  X 10  sec 

D 

2 

Domain  wall  mobility  Is  then  observed  to  follow  quite  closely  a 1/T 

dependence  as  temperature  Is  changed . 

Finally,  using  the  Holsteln-Prlmakof f transformation  techniques,  we 

calculated  relaxation  rates  of  the  uniform  precession  mode  due  to  Intrinsic 

magnon  scattering  processes  In  a uniform  canted  2-sublattlce  YPeOg  system. 

The  calculation  showed  that  the  lowest  order  Intrinsic  processes  that 

lead  to  relaxation  of  the  uniform  mode  are  the  4-magnon  processes,  relax- 

2 

atlon  processes  which  do  have  a T temperature  dependence  In  good  agree- 
ment with  our  experimental  findings.  In  particular,  estimates  of  the 
relaxation  rate  at  room  temperature  also  leads  to  a Bloch  relaxation  time 
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of  T 3 4 X 10"^°  sec 
B 


which  Is  quite  close  to  our  experimental  value, 
considering  the  various  approximations  used  in  this  calculation.  I^istly, 
our  calculations  also  showed  that  it  is  the  scattering  of  2 a-mode  mag- 
nons  into  2 7-mode  magnons , brought  about  by  nondiagonal  terms  of  the 


isotropic  exchange  energy  in  the  YFeO  Hamiltonian  that  constitutes  the 


dominant  contribution  to  the  4-magnon  process  relaxation  rates  that  we 
calculated . 
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ABSTRACT 

A method  Is  described  by  which  GaAs  may  be  encapsulated  to  with- 
stand annealing  temperatures  over  1100“c  using  a double  layered  encap- 
sulant  consisting  of  arsenic  doped  silicon  dioxide  on  top  of  plasma  de-  I 

posited  silicon  nitride.  Samples  encapsulated  in  such  a manner  and 
annealed  show  no  signs  of  mechanical  failure  and  yield  higher  electrical 
activation  of  ion  implanted  selenium  when  compared  with  samples  annealed 
with  silicon  nitride  only.  In  addition,  there  is  no  detectable  outdif- 

fusion  of  Ga  or  As  and  no  detectable  indlffusion  of  Si.  Peak  electrical  i 

19 

activation  of  ion  implanted  Se  has  been  measured  as  1.0  x 10  carriers/ 

3 

cm  for  samples  annealed  at  1100®C.  A first  order  strain  analysis  of 
general  multilayered  systems  is  also  presented,  indicating  possible  im- 
provements on  such  an  encapsulating  procedure  for  GaAs  as  well  as  other 
compound  semiconductors.  The  development  of  this  encapsulant  is  shown 
to  be  essential  for  the  study  of  diffusion  of  ion  Implanted  elements  in 
GaAs  such  as  selenium. 

Electrical  measurements  are  combined  with  the  technique  of  secondary 
ion  mass  spectrometry  (SIMS)  in  order  to  experimentally  anal3'ze  and  cor- 
relate the  diffusion  and  activation  of  ion  implanted  selenium  in  GaAs. 

Pour  chemically  different  species  of  selenium  are  identified:  (1)  sub- 
stitutional selenium,  (2)  interstitial  selenium,  (3)  selenium  complexed 
with  a gallium  vacancy,  and  (4)  precipitated  selenium.  It  is  the  inter- 
actions between  these  four  that  dictates  resulting  redistribution  and 
electrical  activation  of  Implanted  layers.  The  factors  governing  these 
interactions  are  Investigated , and  it  is  demonstrated  that  only  substi- 
tutional selenium  is  a shallow  donor.  In  addition,  it  is  shown  that  the 
species  responsible  for  redistribution  of  impurity  profiles  is  the  sel- 
enium-gallium vacancy  complex.  Precipitated  and  interstitial  selenium 
appear  to  neither  diffuse  nor  act  like  donors  in  GaAs. 

A model  is  developed  which  formalizes  these  observations  in  a set 
of  five  coupled  differential  equations.  By  employing  a minimum  number 
of  simplifying  assumptions,  we  are  able  to  extract  quantitative  predic- 
tions from  this  model  which  accurately  describe  not  only  our  experimental 
results  but  those  of  other  workers. 
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ABSTRACT 

Sensitivity  and  roundoff  errors  can  seriously  limit  the  application 
of  recursive  digital  filters  in  practice,  particularly  when  the  filters 
have  poles  near  z = +1.  A filter  structure,  based  on  digital  incre- 
mental computers,  is  proposed  which  has  low  sensitivity,  good  error 
characteristics,  and  simple  hardware  implementation  for  pole  locations 
close  to  z = +1 . 

The  coefficients  of  the  new  structure  and  the  "conventional"  one 
are  related  by  linear  transformations,  resulting  in  low  sensitivity.  The 
low  roundoff  error  is  a result  of  a closed-loop  scheme  used  in  truncat- 
ing the  outputs  of  the  integrators  and  incremental  multipliers.  The  low  ] 

order  bits  are  not  dropped  but  saved  for  use  at  the  next  epoch.  A sec-  | 

ond-order  digital  filter  Implemented  by  the  new  structure  has  two  arbi-  j 

trary  gain  constants,  which  allow  some  freedom  in  the  design,  and  can  be  | 

used  to  reduce  roundoff  errors  and  limit  cycle  oscillations. 

For  digital  filters  working  in  the  audio-frequency  range,  the  new 
structure  can  be  implemented  without  any  hardware  multiplier.  Also,  for 
a wide  range  of  digital  filters,  excellent  results  can  usually  be  ob- 
tained using  multipliers  with  very  short  word  lengths.  These  can  be 
conveniently  replaced  by  small  size  read-only  memories,  allowing  very 
high  sampling  rates. 

The  accumulation  of  roundoff  errors  is  more  severe  for  digital  fil- 
ters used  for  processing  differentially  pulse-code  modulated  (DPCM)  sig- 
nals. The  new  structure  can  be  directly  used  to  process  DPCM  signals 
with  very  low  roundoff  errors. 

Formulas  for  the  sensitivity  and  roundoff  errors  are  derived  and 
compared  with  those  for  conventional  structures.  Detailed  analysis  of 
limit  cycle  oscillations  is  made  using  a "deterministic"  approach,  and 
the  formulas  obtained  for  limit  cycle  bounds  were  in  perfect  agreement 
with  results  obtained  by  extensive  simulation. 

A design  procedure  is  suggested  to  implement  the  new  filter  struc- 
ture given  the  transfer  function  to  be  realized,  and  the  relations  among 
roundoff  errors,  dynamic  range,  and  the  size  of  the  hardware  multiplier 


1 


ill 


SEL-77-047 


are  discussed.  Implementing  high-order  digital  filters  is  also  consid- 
ered and  simulation  results  are  presented . 

In  short,  the  new  structure  was  found  to  be  "matched"  to  transfer 
functions  which  have  all  their  poles  clustered  around  z = +1. 
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ABSTRACT 


This  paper  concerns  two  Important  Issues  In  the  design  of  optimal 
languages  for  direct  execution  In  an  Interpretive  system:  binding  the 
operand  Identifiers  In  an  executable  Instruction  unit  to  the  arguments 
of  the  routine  Implementing  the  operator  defined  by  that  Instruction; 
and  binding  operand  Identifiers  to  execution  variables.  These  Issues 
are  central  to  the  performance  of  a system,  both  In  space  and  time. 

Historically,  some  form  of  "machine  language"  Is  used  as  the  di- 
rectly executable  medium  for  a computing  system.  These  languages  tradi- 
tionally are  constrained  to  a single  "n-address"  Instruction  format;  this 
leads  to  an  excessive  number  of  "overhead"  Instructions  that  do  nothing 
but  move  values  from  one  storage  resource  to  another  being  Imbedded  In 
the  executable  Instruction  stream.  We  propose  to  reduce  this  overhead 
by  Increasing  the  number  of  Instruction  formats  available  at  the  directly 
executed  language  level. 

Machine  languages  are  also  constricted  with  respect  to  the  manner 
in  which  operands  can  be  "addressed"  within  an  Instruction.  Usually, 
some  form  of  Indexed  base-register  scheme  Is  available,  along  with  a di- 
rect addressing  mechanism  for  a few,  "special"  storage  cells  (l.e.,  reg- 
isters, and  perhpas  the  zeroth  page  of  main  store).  We  propose  a differ- 
ent Identification  mechanism — based  on  the  Contour  Model  of  Johnston. 

Using  our  scheme,  only  N bits  are  needed  to  encode  any  identifier  in  a 
scope  containing  less  than  2^4<N  distinct  identifiers. 

Together,  these  two  results  lead  to  directly  executed  language  de- 
signs which  are  optimal  in  the  sense  that:  (1)  k executable  Instructions 
are  required  to  Implement  a source  statement  containing  k functional 
operators;  (2)  the  space  required  to  represent  the  executable  form  of  a 
source  statement  containing  k distinct  functional  operators  and  v distinct 
variables  approaches  F^k  + N*v — where  there  are  less  than  2**?  distinct 
functional  operators  In  the  scope  of  definition  for  the  source  statement, 
and  less  than  2**N  distinct  variables  In  this  scope;  (3)  the  time  needed 
to  execute  the  representation  of  a source  statement  containing  k functional 
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operators,  d distinct  variables  in  its  domain,  and  r distinct  variables 
in  its  range  approaches  d + r + k;  \diere  time  is  measured  in  memory  ref- 
erences . 
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ABSTRACT 


Intermittent  faults  In  combinational  circuits  may  appear  and  disap- 
pear randomly,  hence  their  detection  requires  many  repeated  applications 
of  test  vectors.  Since  testing  reduces  the  time  available  for  computa- 
tion, it  Is  necessary  to  efficiently  minimize  the  time  required  for  a 
test  while  still  achieving  a high  degree  of  fault  detection. 

This  paper  presents  an  optimal  algorithm  to  detect  intermittent 
failures.  The  algorithm  maximizes  the  probability  of  fault  detection  by 
optimally  choosing  the  input  vector  probabilities. 
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ABSTRACT 


In  present  day  integrated  circuit  technology,  many  faults  cannot  be 
modeled  as  single  faults,  and  the  assumption  of  multiple  faults  is  nec- 
essary. Since  the  class  of  intermittent  failures  constitutes  a major 
part  of  all  field  failures,  a detection  procedure  for  multiple  intermit- 
tent faults  is  extremely  important. 

This  paper  presents  an  optimal  fault  detection  experiment  for  mul- 
tiple intermittent  faults  in  combinational  circuits.  The  optimality 
criterion  used  is  maximizing  the  probability  of  fault  detection.  The 
optimal  experiment  is  achieved  by  applying  input  vectors  with  an  optimal 
probability  distribution. 


ill 


SEL-78-002 


TESTING  FOR  INTERMITTENT 
FAILURES  IN  COMBINATIONAL 
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ABSTRACT 


Since  the  cost  of  detecting  and  Isolating  Intermittent  faults  In 
digital  systems  comprises  a major  part  of  the  total  testing  cost,  effi- 
cient methods  to  detect  them  are  required  . This  paper  proposes  an  op- 
timal algorithm  for  testing  for  intermittent  faults  In  combinational 
circuits  which  are  vital  parts  of  any  computing  system.  The  underlying 
optimality  criterion  is  to  minimize  the  mean  testing  time  and  still 
guarantee  a given  escape  probability.  This  goal  is  achieved  by  apply- 
ing input  vectors  with  an  optimal  probability  distribution. 
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A simulation  technique  has  been  developed  for  the  reliability  eval- 
uation of  arbitrarily  defined  computer  systems.  The  main  simulation 
program  is  written  in  FORTRAN  IV  and  requires  no  changes  to  simulate 
many  different  systems.  The  user  defines  a model  for  a particular  sys- 
tem by  supplying  a set  of  short  FORTRAN  subroutines  and  a specially 
formatted  block  of  numerical  parameters.  The  subroutines  specify  the 
functional  behavior  of  various  subsjrstems  comprising  the  model,  ishile 
the  numerical  parameters  describe  how  the  subsystems  are  interconnected, 
their  time  delays,  what  faults  occur  in  each  one,  etc.  The  main  Simula-  | 

tion  program  uses  this  model  to  perform  a Monte-Carlo-type  evaluation  of 

the  systems'  reliability.  I 

■ 

This  report  supplements  a basic  description  of  the  technique  by  i 

supplying  all  the  details  necessary  for  writing  subroutines,  specifying 
numerical  parameters,  and  using  the  main  simulation  program.  The  simu- 
lation is  event-driven  and  automatically  generates  pseudo-random  faults 
and  time  delays  according  to  parameters  given  by  the  user.  Some  problems 
typical  of  event  simulators,  such  as  ambiguities  arising  from  random 
time-delay  generation,  can  be  solved  by  taking  advantage  of  special  fa- 
cilities built  into  the  simulation  package.  A complete  source  listing 
of  the  main  program  is  included  for  reference. 
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PROGRAMS  WITH  SHARED 
DATA  CLASSES 
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ABSTRACT 


Monitors  are  a valuable  tool  for  organizing  operations  on  shared 
data  in  concurrent  programs.  In  some  cases,  however,  the  mutually  ex- 
clusive procedure  calls  provided  by  monitors  are  overly  restrictive. 

Such  applications  can  be  programmed  using  shared  classes  wiilch  do  not 
enforce  mutual  exclusion.  This  paper  presents  a method  of  verifying 
parallel  programs  containing  shared  classes.  One  first  proves  that  each 
class  procedure  performs  correctly  iidien  executed  by  itself,  then  shows 
that  simultaneous  execution  of  other  class  procedures  can  not  interfere 
with  its  correct  operation.  Once  a class  has  been  verified,  calls  to 
its  procedures  may  be  treated  as  uninterruptible  actions ; this  simpli- 
fies the  proof  of  higher-level  program  components.  Proof  rules  for 
classes  and  procedure  calls  are  given  in  Hoare's  axiomatic  style.  Sev- 
eral examples  are  verified,  including  two  versions  of  the  readers  and 
writers  problem  and  a dynamic  resource  allocator. 
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ABSTRACT 


This  paper  describes  a design  automation  system  for  digital  systems 
that  promotes  a structured  design  approach  both  In  logic  design  phase  and 
In  the  physical  Implementation  phase. 


JSEP  Contract 
N00014-75-C-0601 
Tech . Note  No . 132  , 

November  1977  j 

! 


ABSTRACT 


This  paper  describes  the  initial  design  of  a comprehensive  hierar- 
chical integrated  circuit  design  system,  that  is  currently  being  imple- 
mented at  Stanford  University.  This  system  encourages  the  use  of  struc- 
tured hardware  design  techniques.  It  is  Intended  for  the  design  and 
layout  of  large-scale  integrated  circuits  by  means  of  a combination  of 
manual  and  algorithmic  techniques. 
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ABSTRACT 


This  survey  covers  a wide  range  of  analytical  models  of  dual  redun- 
dant systems.  In  such  systems  two  modules,  each  capable  of  performing 
the  S3rstem  functions,  are  used  to  increase  the  overall  system  reliabil- 
ity. One  module  is  active,  \»hile  the  other  operates  in  a standby  mode. 
The  models  describe  modules  with  and  without  repair  capability.  The 
analysis  of  these  models  gives  us  the  following  information  about  the 
system : 

(1)  the  system  reliability, 

(2)  the  mean  time  before  system  failure,  and 

(3)  the  availability  of  the  system. 

The  models  depend  upon  a variety  of  parameters  used  to  describe  the  sys- 
tem including: 

(1)  the  system  coverage,  i.e.,  the  probability  that  the  sys- 
tem continues  to  function  after  the  active  module  fails, 

(2)  the  probability  distribution  of  the  time  until  recovery 
or  repair  is  completed, 

(3)  the  probability  distribution  of  the  switchover  time, 

(4)  the  failure  rates  of  control  and  switching  mechanisms, 
and 

(5)  the  duration  and  frequency  of  tests  or  checks  of  the 
systems  operation. 

The  analytical  results  from  these  models  are  given  in  this  survey  and 
their  relevancy  to  dual  redundant  computing  systems  is  discussed.  Sev- 
eral computer  systems  are  described  to  demonstrate  the  range  of  applica- 
bility of  dual  redundancy. 
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PROBABILISTIC  TREATMENT  OF 
DIAGNOSIS  IN  DIGITAL  SYSTEMS 
by  M.  L.  Blount 


NSF  Grant 
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Tech.  Note  No.  102 
August  1977 

ABSTRACT 

The  problem  of  self-diagnosis  of  multi-unit  digital  systems  is  con- 
sidered in  this  paper.  Each  unit  can  test  another  unit  and/or  be  tested 
by  another  unit.  The  system  test  consists  of  running  and  compiling  test 
results  from  all  unit  tests  in  the  system.  As  a result  of  the  system 
test,  a decision  is  made  regarding  the  fault  status,  faulty  or  fault-free, 
of  each  unit.  A model  for  performing  automatic  diagnosis  in  an  n-unlt 
system  was  presented  in  the  much  referenced  paper  by  [Preparata,  et  al, 

1967] . 

In  this  paper,  procedures  are  developed  for  calculating  probabilis- 
tic dlagnosability  measures  for  systems  that  can  be  modeled  in  the  graph- 
theoretic  diagnosis  model  of  F>reparata.  The  Preparata  model  is  extended  f 

in  two  important  ways.  Self-testing  units  can  be  modeled;  the  uncertain-  | 

ties  and  imperfections  of  the  testing  process  can  be  modeled.  Algorithms  | 

are  given  for  calculating  1)  the  probability  that  a specified  fault  con-  I 

ditlon  is  correctly  diagnosed  and  2)  the  probability  that  correct  dlag-  | 

nosls  is  performed  when  a system  test  is  done.  Also,  a procedure  is  ^ 

given  for  deriving  a diagnosis  strategy  that  results  in  optimal  system-  ^ 

wide  dlagnosability.  , 
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